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1 Introduction 



Electrodynamics in vacuum or in media with moving boundaries was the subject of numerous studies in the 
XXth century. It is sufhcient to remember that the famous Einstein's paper which gave birth to the relativity 
theory was entitled "Electrodynamics of moving bodies" . The total number of publications in this field 
is enormous, since it includes, in particular, such important problem as radio-location. In this review we 
' confine ourselves to the problem of cavities with (ideal) reflecting moving boundaries. This means that we 

. consider the fields confined in some limited volume, thus leaving aside the problem of the field propagation 

' in the (semi) infinite space and reflection from single boundaries (having in mind that numerous references 

to the publications related to this "full space" problem can be found, e.g., in 1^-^). 

The plan of the chapter is as follows. In the next section we give a brief historical review of the relevant 
studies, both in the fields of classical and quantum electrodynamics, supplying it with an extensive list 
^ i' of the known publications. Although we tried to give a more or less complete list, it is clear that some 

publications have been (undeliberately) missed. To excuse we could mention that for decades different 
^ ' groups of physicists and mathematicians performed the studies in their own fields of interest, not suspecting 

the existence of analogous results found in other areas. We hope that our review will serve to diminish this 
gap essentially. 

It is clear that all results accumulated for several decades cannot be collected in one chapter. Therefore 
in the following sections we have decided to make an emphasis mainly on the detailed exposition of our 
, own results concerning the recently obtained analytical solutions for the cavities with resonantly oscillating 

' boundaries, since we hope that these solutions could be important for the further studies on the problem 

known under the names "Nonstationary Casimir Effect" or "Dynamical Casimir Effect." 

2 Brief history of studies on electrodynamics with moving bound- 
aries 

2.1 Classical fields in cavities with moving boundaries 

The first exact solution of the wave equation {c — I) 



0, (1) 



in a time-dependent domain Q < x < L{t) (where L{t) is the given law of motion of the right boundary), 
satisfying the boundary conditions 

A(0,t)=A(L(t),t) = (2) 
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was obtained by Nicolai ||^ for 

L{t) = Lo{l + at). (3) 

The sohition was interpreted in terms of the transverse vibrations of a string with a variable length. A 
few years later, these results were published in [Q, where the extension to the case of electromagnetic field 
was also made. A similar treatment was given by Havelock ||^ in connection with the problem of radiation 
pressure. Quarter a century later, the one-dimensional wave equation in the time-dependent interval interval 
< X < a + bt was considered in |^ under the name "Spaghetti problem." 

A new wave of interest to the problem of electromagnetic field in a cavity with moving boundaries 
has arisen only in the beginning of 60s, being motivated, partially, by the experiments ]lOt on the "field 
compression," accompanied by the frequency multiplication (for 2.3 times) due to multiple reflections of the 
initial Hqh wave (Ai„ = 10 cm) from the opposite sides of a resonator, one of which was a "plasma piston" 
moving uniformly with the velocity v ^ 2 - lO^cm/s. Kurilko studied linearly polarized electromagnetic 
field between two ideal infinite plates moving with equal constant velocities towards each other. He considered 
consecutive reflections of the waves from each boundary, writing explicit expressions for the finite time-space 
intervals, corresponding to zero, one, two, etc., reflections. Balazs gave a detailed study of the string 
problem with the aid of the method similar to that used by Nicolai and Havelock. Besides considering 

1/2 

the uniformly moving boundary, he has found an exact solution for L{t) = {t^ + l) , and gave some 
graphical method of finding the solution for an arbitrary law of motion L[t). Greenspan JlSf studied the 
one-dimensional string with the uniformly moving right boundary, assuming the boundary condition at the 
left point in the form A(Q,t) = sin(cji). The problem of the "field compression" between two ideal infinite 
moving walls (the same geometry as in p^) was studied by Stetsenko [p^p^. An approximate solution for 
the electromagnetic field in a rectangular waveguide cavity with a uniformly moving boundary was obtained 
in |l^]. In this case one has to solve the equation (c = 1) 

^a-^ - Att = (4) 

with the same boundary condition as in (||). The detailed paper by Baranov and Shirokov can be 
considered, in a sense, as a concluding study of the one-dimensional problem with a uniformly moving 
boundary (although many publications on this subject continued to appear up to last years). The experiments 
on laser cavities with uniformly moving mirrors were described, e.g., in |^"8|-p2|. In these experiments, the 
constant velocity of the mirror varied from 7 cm/s [Q to 400 m/s [ p2[ . 

In short note |2^], Askar'yan has pointed out two possible effects of oscillating surfaces on the elec- 
tromagnetic field inside the (laser) resonator cavities. The first one is the influence of oscillations on the 
generation and intensity of the laser radiation. It was extensively studied in many experiments, devoted, in 
particular, to such problems as the generation of optical pulses |2j,|2^, phase locking of laser modes (where 
the frequencies of the mirror oscillations varied from 50 Hz to 500 KHz |^,^ and 1 MHz ||2^,|3^], see 
the review in pH), or modulation of the laser radiation ]3^,|33|] (in [p3[ the frequencies varied from 17 to 70 
KHz). The theory of these phenomena was considered, e.g., in p3-pq]. 

The second effect predicted by Askar'yan was the field amplification inside the cavity under the parametric 
resonance condition, when the mirror oscillates at twice the field eigenfrequency. It was not observed yet, 
as far as we know, and the main part of the present review is devoted just to the progress in the theoretical 
treatment of this phenomenon achieved in recent years. 

In 1967, Grinberg |^ has proposed a general method of solving the wave equation in the case of an 
arbitrary law of motion of the boundary, based on expanding the solution over the complete set of "instanta- 
neous modes." Krasil'nikov seems to be the first who gave rather detailed study of the electromagnetic 
vibrations in a spherical cavity with an oscillating boundary. The one-dimensional cavity with a resonantly 
oscillating boundary was considered with the aid of the method of characteristics in study , where it was 
shown that the energy is 'pumped' to the high-frequency modes at the expense of the lower-frequency ones. 

A significant contribution was made in a series of papers by Vesnitskii and co-authors. In [^o| he gave 
an exact solution for the problem of a rectangular waveguide with a uniformly moving lateral wall, i.e., for 
the equation and the boundary conditions 

A,, + A,, - Au = 0, Al^^ = Al=L(t) = (5) 
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with L{t) = Lq{1 + at) (in this case, A is the Ey component of the field of the type H„o). A spherical 
resonator, whose radius linearly changed with time, was considered in A general (although implicit, in 
some sense) solution of the one-dimensional problem with an arbitrary law of motion of the boundary was 
given in |^^, where it was shown that a complete set of solutions to the problem (0)-(||) can be expressed 
through the solution of some simple functional equation (see equation (|^) in the next subsection). The 
solutions of the inhomogeneous one-dimensional wave equation for the law of motion L{t) = io(l + oit)'^^ ^ 
with an arbitrary inhomogeneity and arbitrary initial conditions, were given in ||4^ . A family of concrete 
laws of motion admitting simple explicit expressions for the mode functions was found in the framework 
of the inverse problem in p| . The two-dimensional rectangular membrane with a single uniformly moving 
boundary was considered in [ ^Sf . A possibility of the frequency modulation in the waveguide with a slowly 
oscillating boundary was studied in p6| . A review of the results obtained by Vesnitskii and co-authors was 
made in |0 . 

Exact solutions for the waveguide with nonuniformly moving boundary (problem (ph) have been found 



by Barsukov and Grigoryan 1 48 , 49|| . They considered also the electromagnetic resonator with moving bound- 
ary . Similar problems were studied in |5l|, ^ . Periodic solutions of a one-dimensional wave equation 
with homogeneous conditions on moving boundaries were considered in | ]5^ . Transition processes in one- 
dimensional systems with moving boundaries were studied in . Oscillations of a round membrane with a 
uniformly varying radius were considered in ]55| . A nonlinear transformation was applied to solve the inho- 
mogeneous problem of the forced resonance oscillations in a one-dimensional cavity with moving boundaries 
in ||56| . The scaling transformation method, which reduces the problem with moving boundaries to that with 
fixed ones by means of the transformations such as x x/L{t), was considered in |]57|,|5^. 

In the last decade of the century, the problem of vibrating string with moving boundaries was studied 
in Isgj ] (two supports moving toward each other with constant velocities) and |]60|-|6^ (oscillating supports). 
The electromagnetic (better to say, massless scalar) field in one dimensional ideal cavities with periodically 
moving boundaries were considered in ]6^-|69|] . Analytical solutions for the field in circular waveguides with 
(linearly) moving boundaries were obtained recently in [f70|,|7l[|. The "dynamical" field modes in the one- 
dimensional and spherical cavities with uniformly moving boundaries were found once more time in fz^ . The 
same problem for the expanding/contracting ideal spherical cavity, whose radius varies as i?(t) = RovT+at, 
was solved recently in fl^ . (A family of the laws of motion of the boundary, which includes, in particular, 
the dependences such as R{t) = Ro^/T+'at'+'pi^, R{t) ~ Dt + E + F{At + B)^^, and their combinations, 
was considered in the case of the diffusion-type equations in ||7^,[75|], where it was shown that this family 
admits exact solutions of the problem.) 

2.2 Quantum fields in the presence of moving boundaries 

Moore's paper ||76|| seems to be the first one devoted to the problem of quantum fields in cavities mith moving 
boundaries. It was motivated by the studies of a more general problem of the particle (in particular, photon) 
creation in the nonstationary universe |^ and in external intense fields (see, e.g., the books [Q|7^ and 
reviews by Ritus and Nikishov in |^^), which is closely related to the problem of field quantization in the 
spaces with nontrivial (e.g., time-dependent) geometry. Considering a model of the "scalar electrodynamics" 
(when the field depends on a single space coordinate), Moore has found a complete set of solutions to the 
problem (0)-(||) in the form 

An{x, t) — Cn {exp [~iiTnR{t — x)] — exp [~iT:nR{t + x)]} , (6) 

where function i?(^) must satisfy the functional equation 

R{t + L{t)) - R{t ~ L{t)) = 2. (7) 

As a matter of fact, quite similar approach was used in |^-^ for a linear function L{t). Independently, 
equation (^ was obtained by Vesnitskii |^^. Moore's approach was developed in ||7^,|8^,^. However, the 
most efforts were applied to the case of a single mirror (not necessarily plane) moving with a relativistic 
velocity or with a great acceleration, when the effect of particle production becomes significant ]85-10C|. 



The case of two or more boundaries, one of which moves with a const ant r elativistic velocity, was analyzed 



in 101 - 104 (the case of constant relative acceleration was considered in [105|). For the uniform law of motion 
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of the boundary (||) the solution to Eq. (|^), found by many authors cited above, reads (remember that we 
assume c = 1) 

21n|l + , , 

"^"(^^" hrKl + .)/(! 

Evidently, if a — *■ 0, this function goes to i?o(C) = S,/ Lq. For an arbitrary nonrelativistic law of motion one 
can find the solution in the form of the expansion over subsequent time derivatives of the wall displacement. 



We give it in the form obtained in |10£]: 



i?(o = m) - lem + [e - L'io] + • • • , m ^ (q) 

In the special case L{t) = Lo/{l + at), when A(^) = 0, Eq. yields another exact solution, = 
Lq^ + ^a^^). Unfortunately, the expansions such as (^) cannot be used in the long-time limit ^ oo, 
since the terms proportional to the derivatives of A(^) (which are supposed to be small corrections) become 
bigger than the unperturbed term £,X{£,). 

Castagnino and Ferraro have found several solutions of the Moore equation (0) with the aid of the 
inverse method, which was used earlier by Vesnitskii Q]. In this method, one chooses some reasonable 
function i?(^) and determines the corresponding law of motion of the boundary L{t) using the consequence 
of Eq. d), 

.^^^_ R'[t-m]-R'[t + Lit)] 

R'[t-L{t)]+R'[t + L{t)Y ^'"^ 

To solve differential equation (|l0| ) with some simple functions R{£,) is more easy (at least this can be 
done numerically) than to solve the functional equation (|^) for the given function L{t). However, not for 
any simple function i?(^) the dependence L{t) appears admissible from the point of view of physics (the 
velocity may occur greater than the speed of light, or some discontinuities may arise). Actually, the cases 
considered in |]93|| correspond to some monotonous displacements of the mirror from the initial to final 
positions. Typical functions R{£,) used in | p3t were some combinations of £,/Lq and some trigonometric 
functions such as sin(m7rf /Lq). A large list of simple functions R{^) (rational, exponential, logarithmic, 
hyperbolic, trigonometrical and inverse trigonometrical) and corresponding to them functions L{t) was given 
in ]4^ ]. However, no one of these functions can be used in the parametric resonance case. The asymptotical 
solution of the Moore equation in the parametric resonance case L{t) = Lq[1 + esm{TTqt/L())], q = 1, 2, . . ., 
was found in 107| - |llC| ]. For e< ^ 1 it has the form (here Lq — 1) 



R{t) = t-— Ini{ln [1 + C + exp(i7rgi)(l - C)]} , (H) 
C = exp[(-l)''+i7rget] , 

which clearly demonstrates that the asymptotical mode structure in the resonance case is quite different from 



the mode structure inside the cavity with unmoving walls. The solution (11) was improved and generalized to 



the case of two vibrating boundaries in [111, 112|. However, the form of this solution is not very convenient 
for the calculation of various sums giving the mean numbers of photons, energy, etc., so a lot of rather 
sophisticated calculations must be done before the final physical results could be obtained. 

One of the first estimations of the number of photons which could be created from vacuum in a cavity 
whose boundary moves with nonrelativistic velocity have been performed by Rivlin ]113| ] , who considered the 
parametric amplification of the initial vacuum field oscillations in the framework of the classical approach. 
He gave an estimation of the number of created photons M ~ (ewit)^, where e ~ 5L/L is the relative 
amplitude of the variation of the distance between the walls and uji is the fundamental unperturbed field 
eigenfrequency (provided the frequency of the wall vibrations ujw is close to 2cji). A similar estimation 



was given by Sarkar [114|, who used an approximate solution to the Moore equation (|^) in the form of the 



asymptotic series with respect to a small parameter e, found in | 115 | (such solutions were constructed earlier 
in Ip^,^^). However, the tremendous numerical values obtained by Sarkar were quite unrealistic, since he 
used the value of e which was many orders of magnitude higher than those which can be achieved in the 
laboratory under real conditions. Moreover, the oversimplified approach of Rivlin and simple perturbative 
solutions of Sarkar are not valid, as a matter of fact, under the resonance conditions, due to the presence of 
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the secular terms (as it always happens for parametric systems) . Actually, the structure of the 'dynamical' 
modes in the resonance case is completely different (in the most interesting long-time limit) from the simple 
standing waves existing in the cavity with unmoving walls. The same is true for the estimation made by 
Askar'yan |23| in the classical case. He evaluated the average work done by the moving wall on the field as 
A ~ J pvdt, where v = sin(a;„i + </>) is the wall velocity, and p(t) is the radiation pressure. Taking the 
monochromatic dependence p(t) — pa sin^ {2ujit) , he obtained the linear dependence on time in the resonance 
case Lo^ = 2uJi: A ~ povotsincj) + const. However, these evaluations can serve only as some indication that 
in the resonance case the energy of the field can grow at the expense of the mechanical work done by the 
vibrating wall. The real time dependence can be quite different, because the relation p{t) = po sin^(2a;it) 
holds, as a matter of fact, only until ecuit <^ 1, whereas for larger times the effect of the mode reconstruction 
must be taken into account. One of the goals of this review is to demonstrate what happens in reality in the 
resonance case. 

Approximate solutions to the Moore equation, such as (H), w ere used in |ll6| to evaluate corrections to 
the famo us Casimi r attractive force between infinite ideal walls |117| (for the reviews on the Casimir effect 
see, e.g., 118 - 121 1) due to the nonrelativistic motion of the walls. The leading term of these corrections 
turned out proportional to the square of velocity of the boundary (i.e., of the order of {v / cY). The Casimir 
force in the relativistic case was calculated in 0,0|l|,|9O[|l|,||,^|lO|; it depends in the generic case 
not only on the instantaneous velocity, but on the whole time dependence L{t) (through the function i?(^), 
i.e., on the acceleration and other time derivatives). The first calculations of the forces acting on the single 
mirrors moving with nonrelativistic velocities due to the vacuum or thermal fluctuations of the field were 
performed in the framework of the spectral approach (using the fluctuation-dissipation theorem) in |12S] 
(three-dimensional case, the force proportional to the fifth-order derivative of the coordinate) and in |123] 
(one-dimensional model, the force proportional to the third-order derivative of the coordinate). It was 
mentioned that the force could be significantly amplified under the resonance conditions, either in the LC- 
contour [122| or in the Fabry-Perot cavity |123| ] (earlier, such a possibility was discussed in |124]). These 
studies were continued in | 125| -133|, where it was assumed that the velocity of the boundary is perpendicular 
to the surface. The Casimir force between two parallel plates, when their relative velocity is also parallel 
to the surfaces, was considered by Levitov |134|. Later, the theory of "Casimir friction" was developped 
in 1 135 -13!;]. The reviews of these approaches can be found in |14C, 141 1. 

Various quantum effects arising due to the motion of dielectric boundaries, including the modification of 
the Casimir force and creation of photons, both in one and three dimensions, and for different orientations 
of the velocity vector with re spect to the surface, ha ve b een studied in detail in the series of papers by 
Barton and his collaborators |142-148|. In the paper [142| the term Mirror-Induced Radiation (MIR) has 
been introduced. 

Another name, Nonstationary Casimir Effect (NSCE), was introduced earlier in [ 116 | for the class of 
phenomena caused by the reconstr uction of the quantum state of field due to a time dependence of the 
geometrical configuration |14£-152]. Its synonim is the term Dynamical Casimir Effect, which became 
popular after the series of articles by Schwinger |153-157|, in which he tried to explain the phenomenon 
of sonoluminescence by the creation of photons in bubbles with time-dependent radii, oscillating under the 
action of acoustic pressure in the liquids (see a brief discussion of this subject in the last section). 

A possibility of generating the "nonclassical" (in particular, sque ezed) states of the electromagnetic field 
in the cavity with moving walls was pointed out in [106, 114, 124, 15S-16l|. The dynamical Casimir force was 



interpreted as a mechanical signature of the squeezing effect associated with the mirror's motion in [123, 125 



(see also [162 



It was suggested also in [106, 124, 159, 160 [ that a significant amount of photons could be created from 
vacuum even for quite small nonrelativistic velocities of the walls, provided the boundaries of a high-Q cavity 
perform small oscillations at a frequency proportional to some cavity unperturbed eigenfrequency, due to 
an accumulation of small changes in the state of the field for a long time. Indeed, using the asymptotical 
solutions of the Moore equation 



IT|), it was shown in [107-110| that the rate of photon generation in each 
mode becomes constant in the long-time limit, being linearly proportional to the product euji, and that the 
photons are generated in a wide frequency band whose width grows exponentially fast in time |163]. Other 
approximate or exact solutions of the Moore equation for th e specific periodic al time dependences of the 
cavity dimensions were found by different methods in studies [ 111 , 112 , 164- 168 [, which confirmed the effect 
of resonance generation of photons. 
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Moore's approach is based on the decomposition (H) of the field over the mode functions satisfying 
automaticaUy the (one-dimensional) wave equation There exists another approach (proposed in the 
framework of the classical problem as far back as in [37|), when the mode functions are chosen in such a 
way that they satisfy automatically the time-dependent boundary conditions (|^): see Eq. ( |l7| ) below. In 
this case it is possible to describe the behavior of the field with the aid of some effective Hamiltonian, which 
is an infinite-dimensional quadratic form of the boson creation/annihilation operators with time-dependent 
coefficients responsible for the coupling between different modes. Such an approach was considered for the 
first time (in the quantum case) by Razavy and Terning [ 169| -171| (the case of massive field was considered 
in |172|, see also |173]). However, the resulting infinite set of coupled evolution equations for the annihilation 
and creation operators turned out rather complicated in the generic case, and for this reason they were 
treated only perturbatively |I70|. Later, a similar method was used by Calucci [174|, who confined, however. 



only with the case of adiabatically slow motion of the wall, when no photons could be created (recently, this 
adiabatic case was studied in detail in ]I75| ]). An essential progress in the development of the Hamiltonian 
approach (we shall call it also the instantaneous basis method or IBM-method) was achieved after the papers 
by Law |176, 177], who has demonstrated that the effective Hamiltonian can be significantly simplified 
under the resonance conditions. Nonetheless, even the reduced coupled equations of motion resulting from 
the simplified Hamiltonians have been treated for some time either perturbatively (i.e., in the form of 
the Taylor expansion with respect to the time variable) or numerical ly (actua lly, also in the short-time 



limit corr esponding to the initial stage of the process); see the studies [ 176- 181 1 (for the ID cavity model) 
and [182 1 (for a three-dimensional cavity and a waveguide). The general structure of the effective infinite- 
dimensional quadratic Lagrangians and Ha miltonian s arising in the canonical approach to the dynamical 
Casimir effect was a nalyzed and classified in [I83- 185 |. The methods of diagonalization of such Hamiltonians 
were considered in 1 186, 187]. 

The first analytical solutions describing the field inside the ID cavity with resonantly oscillating bound- 
ari es have been found in the simplest cases in il8l - [l90|] . More general solutions have been obtained 
in 1 191 -192]. They hold for any moment of time (provided the amplitude of the wall vibrations is small 
enough; this limitation, however, is quite unessential under realistic conditions). Moreover, these solutions 
enable us not only to calculate the number of photons created from an arbitrary initial state (thus giving, 
for instance, the temperature corrections), but also to account for the effects of detuning from a strict res- 
onance. Besides, they enable to calculate the degree of squeezing in the field quadrature components, to 
find the photon distribution function and the energy density distribution inside the cavity, etc.. Therefore, 
one of the main purposes of this chapter is to give the detailed description of the analytical solutions found 
in 1 188 -192] and to discuss their physical consequences. 

Another goal is to consider the simplest models of the three-dimensional cavity, following the scheme 
given in ]188, 189( ] (for ideal boundaries) and in |194, 195] (for lossy cavities). A more detailed study of the 
quantum properties of the electromagnetic field in rectangular 3D cavities, which takes into account the 
polarization of the field, was performed in ]196,197], but only for the uniform motion of the walls. The 
periodic motion was considered in ]198], also with account of the polarization and the influence of all three 
dimensions, but in the framework of some approximations equivalent to the short-time limit. The case of 
a three-dimensional rectangular cavity divided in two parts by an ideal mirror, which suddenly disappears, 
was considered in ]19£]. 

The name M otion Induced Radiation was given to the effect of ra diati on emission outside the cavity with 
vibrating walls ]200|. Using the spectral approach, the authors of ]200] showed that the radiation can be 
essentially enhanced under the resonance conditions (by the orders of magnitude, comparing with the case 
of a single mirror). The problem of the photon generation by a single perfectly reflecting mirror performing 
a bounded nonrelativistic motion was studied in ] 201 , ^02[ , where the effects of polarization have been taken 



into account, and the spectral and angular distributions of the emitted photons have been found. Arbitrary 
space-time deformations of a single moving mirror have been treated in ]203, 204] with the aid of the path- 
integral approach. The same approach was used in the case of a cavity with deformable perfectly reflecting 
boundaries in ]205]. 

The creation of photons or specific (e.g., squeezed) states of the electromagnetic field due to the motion 
of some effective mirrors made of the free electrons moving with (ultra)relativistic velocities was considered 
in |20(;,207]. Another kind of "effective moving mirrors" consisting of the electron-hole plasma generated in 
semiconductors under the action of powerful laser pulses was suggested in [20g,209|. 
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The influence of te mper ature on the dynamical Casimir effect was evaluated inj 110 , 210|| ; a more detailed 
analysis was given in 211 ( in the fr amework of the Hamiltonian approach), ||212|| (using the thermofield 
dynamics), and especially in |192, 193 1. The energy density distribution inside the cavity under the resonance 
conditions is no more uniform, on the contrary, the main part of the energy is concentrated in several sharp 
peaks which move from one boundary to another, becoming more and more narrow with the course of time. 
This effect was studied in the framework of numerical calculations in [112, 164- 164 168| (earlier, it was 
discussed in for the classical field). The analytical form of the pulses, including their "fine structure" 
in the case of initial states different from the vacuum or thermal ones, was found in [193|. A similar pulse 
structure of radiation emitted from the high finesse vibrating cavity with partially transparent mirrors was 
studied in plO| , p3| . 

The evolution of classical fields in the cavities filled in with media whose dielectric properties vary in 
time was considered, e.g., as far back as in [214|. Yablonovitch [215| proposed to use a medium with a 
rapidly decreasing in time refractive index ("plasma window") to simulate the so-called Unruh effect 216|, 
i.e., creation of quanta in an accelerated frame of reference. More rigorous and detailed studies of quantum 
phenomena in nonstationary (deformed) media have been performed in |15£, 160, ^17 - p29 |. The case when 
the dielectric constant changes simultaneously with the distance between mirrors (in one dimension) was 
considered in [230, 231 [. A comparison of the spectra of photons created due to the motion of mirrors and 
due to the time variations of the dielectric permeability was made in ]232[ . An analog of the nonstationary 
Casimir effect in the supe rflui d ^He, namely, the friction force on the moving interface between two different 
phases, was discussed in [233[. 

The problem of the interaction between the electromagnetic field created due to the NSCE and various 
detectors (harmonic oscillators, two-level systems, R ydberg atoms , etc .,) placed inside the cavity with moving 
walls was studied by means of different methods in 18^ , 189, ^34 -23^1. It is also discussed in this chapter. 

The first experiments on the interaction between the powerful laser radiation and freely suspended light 
mirrors have been reported in [239, 240|[ , where the effect of the optical bistability, similar to that observed 
usually in the so-called Kerr media, was observed. The first theoretical study of this phenomenon in the 
cavities whose walls can move under the action of the radiation pressure force was given in |241[. This 
subject received much attention for the last decade in connection with different problems, such as, e.g., 
the attenuation or elimination of noise in interferometers p4^|25^ . It is important, in particular, for the 
gravitational wave detectors and for the general problem of measuring weak forces acting on a quantum 
system: see, e.g., [254, 255| [ for more details. Another possible application could be the generation of the 
so-called "nonclassical states" of the field and the mirror itself (when it is considered as a quantum object, 
too) [ ^53 , 256 -258 1 . The stability of such states, in turn, is closely related to the general problem of decoher- 
ence in quantum mechanics, and one of the mechanisms which can destroy the coherence of pure quantum 
superpositions is just the dynamical Casimir effect [259, 260 [. The back reaction of the dynamical Casimir 
effect was recently studied in [261[. The influence of fluctuations of the positions of the walls on the field 
inside the cavity was considered in [115, 205, 262 [, whereas the Brownian motion of the walls due to the field 
fluctuations was studied in [|140,^63 [. The role of the dynamical Casimir effect in the cosmological problems 
was studied recently in [ 264 , 26£ ] . 



3 ID cavity with oscillating boundaries 

Let us start with the case of a single space dimension. Consider a cavity formed by two infinite ideal plates 
moving in accordance with the prescribed laws 

Xleftii) = U{t), Xright{t) = u(t) -f- L{t) 

where L{t) > is the time dependent length of the cavity. Taking into account only the electromagnetic 
modes whose vector potential is directed along z-axis ("scalar electrodynamics" [|7^]), one can write down 
the field operator in the Heisenberg representation A{x, t) at t < (when both the plates were at rest at the 
positions xieft — and Xright — Lq) as (we assume c — h = 1) 

oo 

A ^ — A i . TITTX f" / \ 1 /-I r\\ 

Ain = 2} ^sm— — 0„exp (-iWnt) + h.c. (12) 
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where &„ means the usual annihilation photon operator and LOn = im/LQ. The choice of coefhcients in 
equation (|l2|) corresponds to the standard form of the field Hamiltonian 



H 



1 



dx 



OA 
'dt 



OA 
dx 



°° „ 1 



For t > the field operator can be written as 



h.c. 



(13) 



(14) 



To find the explicit form of functions ^^'^''{x,t), n — 1,2, . . ., one should take into account that the field 
operator must satisfy 

i) the wave equation (Q), 

ii) the boundary conditions (H) or their generalization 



A{u{t), t) = A{u{t) + L{t), t) = 0, 
iii) the initial condition (p^), which is equivalent to 

^/i*"^ (x, t < 0) = sin exp (-iw„t) 



(15) 



(16) 



Following the approach of Refs. 170, 174, 176, 191| we expand the function '4>^''^\x,t) in a series with 
respect to the instantaneous basis: 



^^-Hx,t > 0) = Egi"^(^)^^sin f ^M^^^ 

with the initial conditions 



fc=i 



Lit) 



n= 1,2,. 



(17) 



Qi"^(0) = 4„, Oi"^(0) = -zc^„4„, fc,n=l,2,... 

This way we satisfy automatically both the boundary conditions (^) and the initial condition ( p^ . Putting 
expression (17) into the wave equation (Q), one can arrive after some algebra at an infinite set of coupled 
differential equations p80|,p4|p9l| 



9%) 



where 



LUk{t) = kTr/L{t) 

and the time dependent antisymmetric coefficients gkj (t) read (for j ^ k) 



9kj 



-9jk = (-1)" 



2k j (L + uckj ) 



if -k^) Lit) 



(18) 



(19) 



For u — (the left wall at rest) the equations like (p^-(pj|) were derived in |174, 177|. 

If the wall comes back to its initial position Lq after some interval of time T, then the right-hand side of 
equation (nSh disappears, so at t > T one can write 



(20) 
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where f^"'' and 77^"' are some constant complex coefBcicnts. Consequently, at t > T the initial annihilation 
operators &„ cease to be "physical" , due to the contribution of the terms with "incorrect signs" in the expo- 
nentials exp{iujkt). Introducing a new set of "physical" operators dm and a]„, which give the decomposition 
of the vector potential operator at t > T in the form analogous to (E2|), 



00 2 



(21) 



V — ^ ^ 

A{x, t) — — -= sin (irnx/Lo) 

n— 1 ^ 

one can easily check that the two sets of operators are related by means of the Bogoliubov transformation 

00 

am = X! (^n*^"™ + ^if^nm) ' m = 1, 2, . . . (22) 

n=l 

with the coefficients 



«n™ = A/-4"\ /3nm = J-r?'?). (23) 

V n y n 

The unitarity of the transformation ( |2^ ) implies the following constraints: 

oc oo 



n 

m— 1 



oo oo 

E («:™««. - /9:m/3n.) = E ^ (eL"^*ej"^ - ^/i?)*.?? ') = '5™, (25) 

n—l n—1 
00 00 ^ 

E ^f^nm^n, " /3:.an™) = E " {^^n*^"^ " ^i"^*^"^) = (26) 



n 

n—l n—l 



The mean number of photons in the mth mode equals the average value of the operator d]„dm in the 
initial state |in) (remember that we use the Heisenberg picture), since just this operator has a physical 
meaning at i > T: 

J\f„i = (in|aj,ja,„|in) 

= Ei/5"'"i' + E[( 

+ /3:™/3fem)(^\fe/c)+2Re {bnbk] 

n n,k 

oo oc 



1 , , vnfc 

n=l n,fc=l 
oo 

+2Re E -^^t\\nKKhK). (27) 

The first sum in the right-hand sides of the relations above describes the effect of the photon creation from 
vacuum due to the NSCE, while the other sums are different from zero only in the case of a nonvacuum 
initial state of the field. 

To find the coefficients and 77^"^ one has to solve an infinite set of coupled equations (^ (fc = 1, 2, . . .) 
with time-dependent coefficients, moreover, each equation also contains an infinite number of terms. However, 
the problem can be essentially simplified, if the walls perform small oscillations at the frequency cuio close to 
some unperturbed field eigenfrequency: 

L{t) = Lo (1 + £l sin [puji{l + S)t]) , u{t) = EuLq sin [pwi(l + S)t + ip] . 

Assuming |£l|, |£u| ~ e <C 1, it is natural to look for the solutions of equation in the form similar to 
(§), 

Q<^'\t) = ^(")e— (1+^)* + ,7(")e*"'=(i+^'*, (28) 
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but now we allow the coefficients and 77^"'' to be slowly varying functions of time. The further procedure 
is well known in the theory of parametrically excited systems |266- ^68 |. First we put expression ( p8[ ) into 
equation ([l^) and neglect the terms |, i] (having in mind that S,,t] ^ s, while 1^, 77 ~ e^), as well as the terms 
proportional to ~ ~ e^. Multiplying the resulting equation for Qk by the factors exp[ia;fc(l + 
and exp [— zcjfe(l + S)t] and performing averaging over fast oscillations with the frequencies proportional to 
LUk (since the functions ^,77 practically do not change their values at the time scale of 27r/cjfc) one can verify 
that only the terms with the difference j — k = ±p survive in the right-hand side. Consequently, for even 
values of p the term u in gkj{t) does not make any contribution to the simplified equations of motion, thus 
only the rate of change of the cavity length L/ is important in this case. On the contrary, if p is an odd 
number, then the field evolution depends on the velocity of the centre of the cavity Vc — u + L/2 and does 
not depen d on L alone. These interference effects were discussed in the short time limit etoit <C 1 in |l8Ct] 
(see also |200]). We assume hereafter that u — Q (i.e. that the left wall is at rest), since this assumption 
does not change anything if p is an even number, whereas one should simply replace L/Lq by 2vc/Lo if p is 
an odd number. 

The final equations for the coefficients ^^"■^ and 77^"'' contain only three terms with simple time independent 
coefficients in the right-hand sides: 

_d_ („) 
dr^'^ 



(n) 



i-ir 
{-ir 



{k- 
[k- 



(n) 



2ijkri'j^\ 



(29) 
(30) 



The dimensionless parameters r (a "slow" time) and 7 read (e = el) 



The initial conditions are 



ei"Ho) 



Skr, 



7 = S/e. 



0. 



(31) 



(32) 



Note, however, that uncoupled equations (^9[)-(^0|) hold only for k > p. This means that they describe the 
evolution of all the Bogoliubov coefficients only if p = 1. Then all the functions 'n'k\t) are identically equal 
to zero due to the initial conditions (p2[), consequently, no photon can be created from vacuum. Moreover, 
in the next section we show that the total number of photons (but not the total energy) is an integral of 
motion in this specific case. 



4 "Semi-resonance" case (p = 1) 

If j3 = 1, one has to solve the set of equations {k, n — 1,2, . . .) 

= {k- l)eS - (fc + l)ea + 277fc^r^ . (33) 
An immediate consequence of these equations and the condition (0) = Skn is the identity 

J2me„:\r)^i':>{r)^nS,,k, (34) 

rn 

which is nothing but the unitarity condition of the Bogoliubov transformation in this special case. Taking 
into account this identity, one can easily verify that the total average number of photons in all modes is 
conserved in time: 

J\f=J2 -^d"^*d'^(in|6t 6fe|in) = ^(in|6t &„|in). (35) 

nkm n 

A similar phenomenon in the classical case was discussed in [^4| , whereas the quantum case was considered 
in 176 and especially in |19C]. Also, using Eqs. ( |33| ) and ( ^4| ) one can verify that the total energy (normalized 
by uji) 

nkm 
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satisfies the simple equation (hereafter overdots mean the differentiation with respect to the dimensionless 
time r) 

E = ^c?E + 4725(0) - 27lm(gi), (36) 
where 

(37) 



(38) 



The quantum averaging is performed over the initial state of the field (no matter pure or mixed). The initial 
values of the total energy and its first derivative (with respect to r) are given by 



f(0) = ^n(6tfe„), f(0) = Re(^i). 



(39) 



Consequently, the solution to equation (|3^) can be expressed as 

2 sinh^ (ar 



^ 7 , „ ^^ ^ ,„ , sinh(2ar) 
£(0) - ^Im(gi)J + Re(gi) 



(40) 



One can easily prove that |£(0)| < f (0). Thus the total energy grows exponentially when t ^ 1 (provided 
7 < 1), although it can decrease at r ^ 1, if £(0) < 0. Since the total number of photons is constant, such 
a behaviour is explained by the effect of pumping the highest modes at the expense of the lowest ones (in 
the classical case this effect was noticed in Q). 

We see that the total energy can be found without any knowledge of the Bogoliubov coefficients. However, 
these coefficients are necessary, if one wants to know the distribution of the energy or the mean photon 
numbers over the modes. To solve the infinite set of equations (p3[) we introduce the generating junction 



(41) 



k=l 



where z is an auxiliary variable. Using the relation kz'^ — z(dz'^/dz) one obtains the first-order partial 
differential equation 



(z^ - 1 + 2i7z) 



dr dz 
whose solution satisfying the initial condition X'^")(0, z) — z" reads 

zg{T) - S{t) 



X(")(z,r) 



5*(t) - zS{t) 



+ I d"^(:c)dx 



where 

S{t) = sinh(aT)/a , ^(t) = cosh(aT) + ijS{T) . 
Differentiating (^) over z we find 



nhSir)] 



n-l 



(42) 

(43) 
(44) 
(45) 



[5*(r)]"+i 

Putting this expression into the integral in the right-hand side of equation (^3|) we arrive at the final form 
of the generating function 

X(")(z,t) 



r) = 


' zgir) - S{t)' 


n 


\-s{T^ 




[g*{T)-zSiT)\ 







(46) 

which satisfies automatically the necessary boundary condition X^'"-\t,0) = 0. The right-hand side of ( ^ ) 
can be expanded into the power series of z with the aid of the formula ( ]269( |, vol. 3, section 19.6, equation 
(16)) 



(i-ty-'ii-t + xty 



A — ^ 777,1 



m— 
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where F(a,b] c] x) is the Gauss hypergeometric function, and (c)^ = r(c + k)/T{c). In turn, the function 
{c)mF{—m,b;c;x) with an integer m is reduced to the Jacobi polynomial in accordance with the formula 
( ||269| , vol. 2, section 10.8, equation (16)) 

(c)™F(-m, 6; c; x) = m!(-l)"Pi''-'"-^- (2a; - 1). 

Consequently, 



{l-t)''-''{l-t + xt)-^ ^ ^(-t)™Pi^-"-'=''=-i)(2a;- 1) (47) 



m=0 



and the coefficient -* (r) reads 

' (^) = (-k)""'"A"+"p/„"-'"' (l - 2^2) (48) 

where 

^ N S(t) 
K r = ^= = ^ ^ ^ 49) 



Mr) = V5(T)/ff*(T) = v/l-7'^'+i7«, |A| = 1. (50) 

The form (|4^) is useful for n > m. To find a convenient formula in the case of n < w we introduce the 
two-dimensional generating fmiction 

oo oo oo 

X{r,z,y) = ^ J^z-y-d-Hr) = 5]x(")(z,r)y" 

m—1 n—1 n—1 

(51) 



[5* (r) + 2/5(r)] [g* (r) - g(r)y^ + S{t) [y - z)Y 
The coefficient at z™ in (pi]) yields another one-dimensional generating function 



X fr u)- Vu"f(")M-wi^^^^^i^i^ll^ (52) 
^m[T, y) - 2_^y (^rj - y ^^^^^^ ^ yS'(r)]™+i ' ^ ' 

Then equation (^^ results in the expression 

e^^) = (1 - 1) _ 2«;2) . (53) 

Note that the functions S{t), cosh(ar) and k{t) are real for any value of 7. For 7 > 1 it is convenient to 
use instead of (E3) the equivalent expressions in terms of the trigonometrical functions: 



S{t) — sin(aT)/a , ^(t) = cos(aT) + i"fS{T), a — \J^'^ — \ (54) 
In the special case 7=1 one has S{t) — r and gir) = 1 + ir. In particular. 



The knowledge of the two-dimensional generating function enables to verify the unitarity condition (|2^) . 
Consider the product X*(r, zi, j/i)X(t, Z2, 2/2), which is a four-variable generating function for the products 
S,m^*£,i''^ ■ Taking yi — ^/uexp(iip), j/j = y/ueyi'p{—iip) and integrating over ip from to 27r one obtains a 
three-variable generating function z™z!2u"^m'^*^i^\ Dividing it by u and integrating the ratio over u 
from to 1 one arrives finally at the relation 

00 00 -. 

E ^r4^e)*^l"^=-ln(l-;.^2) = E^(^^2)\ (56) 

n,m,l — l k—1 

which is equivalent to the special case of ( p5| ) for 7]^'' ^ 0: 
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4.1 Examples 



Suppose that initially there was a single excited mode labeled with an index n. Due to the linearity of 
the process one may assume that the mean number of photons in this mode was = 1. Then the mean 
occupation number of the m-th mode at t > equals 



(n) 



m 
n 



(1 - k.^)k'^ 



-11 p(ni-n, 1) 



1 - 



(58) 



where k is given by ([49|). For example, in the special case 7 = we have 

\2m-2 



« r(i) ^ mftanhr)^ 
(coshr)4 



2(coshr)^ 



1) - (to + 1) tanh^ r] 



The maximum of function A/'m (r) is achieved at sinh Tmax = \J (jn — l)/2. For m 3> 1 it equals M)ri ijraax) 



(1)/ 



4/ (me**). For a fixed value of r 3> 1, the occupation number distribution Mm reaches its maximum at 
.j^ = cosh^ r, and N^l. 



m 



2 \ ^1 

e cosher < 1. The maximum of the energy distribution is shifted to 



,(1) 



2 • • • • fi") 

2 cosh r, and its value is not decreased with time: £max 



4/e^. This explains the 



the right, TO^^^g 
exponential growth of the total energy. 

Although formula (Kq) seems asymmetrical with respect to the indices to and n, actually the relation 



holds. To prove it we calculate the generatmg function 

00 



(59) 



(60) 



m,n— 1 



It is related to the function X{z,y) ( pl| ) as follows 



Q[u,v) = V — 
av 



dr 



2n /.2ir 



Having performed all the calculations we arrive at the expression 

1 + UV ~ {u -\- v) 



2Q{u,v) 



I [1 + — k2(u + v)]^ — 4uw(l — K^)^! 



1/2 



- 1. 



(61) 



Then ( p9|) is a consequence of the relation Q(m, v) = (5(w, w). 

The initial stage of the evolution of Mm\T) does not depend on the detuning parameter 7, since the 
principal term of the expansion of (|5q) with respect to r yields 



AAi^,(r->0) 



n± q 



n{n ± 1) . . . (n ± q =F 1) 



r29 



However, the further evolution is sensitive to the value of 7. If 7 < 1, then the function Afm^r) has many 
maxima and minima (especially for large values of to and n), but finally it decreases asymptotically as 
TOna^/ cosh'*(aT). On the contrary, if 7 > 1, then the function A/'i?'' (r) is periodic with the period 7r/a, and 
it turns into zero for r = kn/d, k = 1,2,... (excepting the case m — n). The magnitude of the coefficient 
■M'm\T) decreases approximately as 7^2!™-"! Jqj. ^ 1 

Now let us assume for simplicity that 7 = 0. Then Eq. ( ^3|) can be represented in the equivalent form 



in) _ 



(tanhr)"- 
cosh^ T 



(-1)"-^F l-n,TO+l;2; 



1 



cosh r 



(62) 
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If m ^ n ~ C(l), then {m)/, « m'^ {k < n), so the Gauss hypergeometric function in Eq. ( |62|) can be 
replaced by the confluent hypergeometric function with a negative integral first index, which is reduced 
to the associated Laguerre polynomial p69[ of the scaled variable /i = m/ cosher. Using the 

approximation (tanhr)^™ « exp {—m/ cosh r) valid for r ^ 1, we arrive at a simplified expression 

- ™AClr) = -M^e-^ ' , (63) 

n L J 

describing the energy distribution over the modes with large number to. 

The fluctuations of the occupation numbers can be calculated with the aid of the formula (again for 

7-0) 

n'^ L 1 L J 

which is an immediate consequence of Eq. (|57|). Consequently, the type of the photon statistics (sub- or 
super-Poissonian) is conserved. In particular, if the initial mode was in a coherent state, then all other 
modes will be excited in the coherent states, too. 

The assumption on the equidistant eigenmode spectrum can be justifled to certain extent for the longitu- 
dinal modes of a Fabry-Perot resonator with perfect mirrors, if the order of interference (the mode number) 
is high enough. So let us suppose that initially some single mode with n 1 was excited. Although the 
lowest modes are not equidistant in this case, for a limited period of time in the beginning of the process 
(until T <C log(2n)), the evolution of the chosen mode and its neighbours can be described by the solutions 
obtained above, since the influence of the "remote" modes becomes essential at sufficiently large times. Then 
the populations of the modes exhibit strong oscillations, since they are proportional to the squares of the 
Jacobi polynomials of large degrees. The asymptotics of these polynomials ( p69[ , Eq. (10.14.10)) yields 

■^1"^ ~ -7 7 cos^ ((to + n)ip - [2|to - rtl + 1]-| , sin^s = tanhr. 

n7rsm(2(y9) L 4 J 

This formula holds provided that n ^ |to — n| ~ 0(1) and nsin(2(p) ^ 1. For r ^ 1 we observe fast 
oscillations, which amplitude is modulated with the period tt/t, whereas for r ^ 1 we have more slow 

periodic variations of Mm"^ as the function of to with the period ire^ /2. 

An interesting problem is the field evolution in a cavity which was initially in the equilibrium state 
at a finite temperature, when the initial occupation numbers were given by the Planck distribution Un = 
[exp(/3n) — 1]^^. Let us consider two limit cases. The first one corresponds to the low-temperature approx- 
imation Vn = eKp(—Pn). Then the occupation number of the TO-th mode is nothing but the coefficient at 
V™ in the expansion (p3l|) with u = exp(— /3). Using the well known generating function of the Legendre 
polynomials Pmiz) ( [ ^69[ , Eq. (10.10.39)), one can obtain the following expression (for 7 = 0): 

^"^^ 2 ( r-e-v ) iPrn{x)+P„,-i{x)], p = tanh(r). 



e 



(e-'3-p2)(l_g-/3^2) ■ 

In particular, A/'/'^^ = e~''(coshr)~'' (l — tanh^ r) ^ . 

In the special case of a cavity filled in with a high-temperature thermal radiation, the initial distribu- 
tion over modes reads 1^,1(0) = Q/n, constant O being proportional to the temperature. Then Mni^^ = 
'Ylin^n{®)J^rn^ ■ This sum is nothing but multiplied by the coefficient at in the Taylor expansion of 
the function 



Q{v) = / — Q(u,w) =ln- ^ 

Jo u 1-v 

Thus we have 

f,if}-TOAA,{f> = e(l-Kr)p™). 

We see that the resonance vibrations of the wall cause an effective cooling of the lowest electromagnetic 
modes (provided I7I < 1). The total number of quanta and the total energy in this example are formally 
infinite, due to the equipartition law of the classical statistical mechanics. In reality both these quantities 
are finite, since i^„(0) < Q/n at n 00 due to the quantum corrections. 
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5 Generic resonance case p > 2 

If p > 2, we have p ~ 1 pair of coupled equations for the coefficients with lower indices 1 < fc < p — 1 



in) 



in) 



dr 



(n) 



i-ir \ik+p)viZ--ip-k)i% 



(64) 
(65) 



In tliis case some functions T]k^\t) iiot equal to zero at t > 0, thus we have the effect of photon creation 
from the vacuum. 

It is convenient to introduce a new set of coefficients p^"'' , whose lower indices run over all integers from 
— oo to oo: 

0, k^O (66) 



(") 
Pk = 



^k ' 



A: < 



Then one can verify that equations (p9|)-(|30|) and (|64[)-(|65[) can be combined in a single set of equation 
(fc = ±l,±2,...) PI 



--pi") = (-i)^ 



(") 



An) 



dr 



with the initial conditions {n — 1,2,. 



2z7fcpi") 



(67) 



•) 



pi"^(0)-4n. (68) 

A remarkable feature of the set of equations ( |67|) is that its solutions satisfy exactly the unitarity conditions 
(p^)-(p6|) (although the coefficients and 77^"'' introduced via equation (^8|) have additional phase factors in 
comparison with the coefficients defined in equation (|2^), these phases do not affect the identities concerned), 
which can be rewritten as 



mpl^'>* pI^^ = nSnk , n,k^ 1,2,... 



OO 



(n)* (n) 
P-,n P-j 



= Srnj , m,j = l,2,... 



(69) 
(70) 



= 0, m,j=:l,2,. 



(71) 



For example, calculating the derivative / = (d/dr) J2m=-oo "^P^m* p^m with the aid of equation (67) and its 
complex conjugated counterpart one can easily verify that / = 0. Then the value of the right-hand side of 



( |69| ) is a consequence of the initial conditions ( |68| ) . The identities ([70|) and (|7l|) can be verified in a similar 
way, if one uses instead of (^) the recurrence relations between the coefficients p^m with the same lower 
index m but with different upper indices: see equatio ns (p7| ) and (|88|). 

Due to the initial conditions (68) the solutions to (|67[) satisfy the relation 



(72) 



j,/c-0,l,...,p-l, m = 0,±l,±2, 

(n) 



n = 0,1,2, 



Consequently, the nonzero coefficients pm form p independent subsets 

= (j+1P) 
^k ^ j+kp 



(73) 



i-0,l,...,p-l. 



0,1,2, 



A: = 0,±l,±2,. 
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The subset 2/^''°'' is distinguished, because yl?'°^ = for A; < and the upper index q begins at g = 1. 
One can verify that the functions ym'^\T) with m > 1 are given by the formulae for ^^'(t) found in the 
preceding section, provided one replaces t by (— l)^pT and 7 by (—1)^7, whereas ym'^\T) = for to < 0. 

In the generic case j 7^ it is reasonable to introduce a generating function in the form of the Laurent 
series of an auxiliary variable z 



m— — 00 



(74) 



since the lower index of the coefficient yin runs over all integers from —00 to 00. One can verify that the 
function (^3) satisfies the homogeneous equation 



dr 



1 



cr I z 



2i7 



j+pz— <J={-1)P. 



(75) 



The solution to (|75|) satisfying the initial condition i?("'^'(z,0) — z" reads 



i?("'j')(2 



-j/p 



zg{pT) + aS{pT) 
g*{pT) + z(tS(jpt) 



n+j/p 



(76) 



where the functions S{t) and ^(r) were defined in (44). The coefhcients of the Laurent series ( [74[ ) can be 
calculated with the aid of the Cauchy formula 



2m 



Z'' 



(77) 



where the closed curve C rounds the point z = in the complex plane in the counterclockwise direction. 
Making a scale transformation one can reduce the integral (^^ with the integrand (76) to the integral 
representation of the Gauss hypergeometric function ( ]269|], vol 1, section 2.1.3) 



(78) 



-»r(c) cxp(-^7^b) t^-\l - tf-^-^ 

^' ^ 2sin(.fe)r(c-6)r(6) (1 - t.Y 



where Re(c — 6) > 0, 6 7^ 1, 2, 3, . . ., and the integration contour begins at the point t = 1 and passes around 
the point t = in the positive direction. After some algebra one can obtain the expression 



_ ^ r (-TO - j/p) r (1 + n + j/p) sin [tt (to + j/p)] 
ttF (1 + n — to) 



n — m: K 



(79) 



We assume hereafter k = k{pt) and A — X{pt), the functions k{x) and X{x) being defined as in (|49| ) and 
(^). Using the known formula 

r(-z)sin(7rz) = -7r/r(z + l) (80) 
one can eliminate the gamma- function of a negative argument: 

(„^j) _ r (1 + Tl + j/p) (cr^)"-'»A"+"+^j/P 



r (1 + TO + j/p) r (1 + n — to) 

X F (n + j/p , —TO — j/p ; 1 + n — TO ; k^) 



(81) 



The form (81) gives an explicit expression for the coefficient with < to < n. Moreover, it clearly 

shows the fulfilment of the initial condition ym''^\T = 0) = <^mn- Transforming the hypergeometric function 
with the aid of the formula |26!; , 270 



lim 

c — *—n 



F{a,b;c;x) {a)n+i{b)n+ix 



n+l 



r(c) 



(n + l)! 



-F{a + n + l,b + n + 1; n + 2; x) 
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{n = 0, 1, 2, . . .) and the identity (|8C|) one obtains an equivalent expression 

(„^^-) _ r (to + j/p) (-crK)'"-"A"+"+^J/P 



T (n + j/p) r (1 + TO — n) 
F (to + j/p , ~n — j/p ; 1 + TO — n ; k^) 



(82) 



which gives a convenient form of the coefficient S.j'^j^^ for m> n. Formula (^9|) with negative values of the 
lower index gives an explicit expression for the nonzero coefficients J^p^"^"*' (fc > 1, n > 0): 

ipn+j) ^ _ r (fc - j/p) r (1 + n + j/p) sin [tt (fc - j/p)] 



■7TT{l + n + k) 

X (ctk)"+'=A"-'=+^^'/pF {n + j/p,k^j/p-,l + n + k-, k^) 



(83) 



Note that the expressions (]8l|)-(|83|) are valid for j — 0, too. In this case they coincide with the formulae 
obtained in the preceding section. Formulas (^)-(|83|) immediately give the short-time behaviour of the 
Bogoliubov coefficients at r — )■ 0: it is sufficient to put k « pr, A « 1 and to replace the hypergeometric 
functions by 1. In this limit the detuning parameter 7 drops out of the expressions (in the leading terms of 
the Taylor expansions). 

At T — !■ 00 we have the following asymptotics of the functions k{pt) and X{pt) (if 7 < 1) 



1 



l-^S-'{pr)^l, A 



Then equation (£9 



at a; < 1 |269, 270| 



a + tj, T — !■ 00. 

together with the known asymptotics of the hypergeometric function F{a, b; a+b+l; 1—x) 



F{a,b;a + b+l;l- x) 



r(a + 6+ 1) 



r(a + i)r(6 + i) 

lead to the asymptotical expression for the Bogoliubov coefficients 



[1 + abxlnix) + 0{x)] 



(84) 



(r»l) = ""["("^+.^'/f^] (a + Z7r 



7:{m + j/p) 
/"mn 



i-\-n-\-2j/p^7i—r 



i + o{ 
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For 7 < 1 the correction has an order TOnTexp(— 2apr), while for 7 = 1 it has an order mriln(T)/T2. 
One can verify that the generating function (|7|) satisfies the recurrence relation 



(85) 



dr 



{j + qp) {<7 R^'^'^'^^ - + 2i7i?(«'^)} 



(86) 



Its immediate consequence is an analogous relation for the Bogoliubov coefficients with the same lower 
indices: 

d 



dr 



rm 



r7n rm 



(87) 



Equation ( |8^ ) is valid for n > p (when q > 1 and j > 1 in (p6|)), since the coefficients pm"^ are not defined 
when n < 0. However, using the chain of identities 



j/p 


' S + gz' 


-I 






-9*/z] 




[g* + Sz\ 






. 9- 


VS/z 



i z ^-^ 



k— — oo 



7 - E 



{0,p~j)*^k 



one can obtain the first p — I recurrence relations 
d 



dr 



Pn 



(p-n)* (p+n) 
r—rn rm. 



n = 1,2,. 



,p-l. 
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To treat the special case n = p {it corresponds to the distinguished subset with j = 0) one should take into 
account that = 1, which means formally that = Smo- So the last recurrence relation reads 

= p { -apt^P) + 2i7p(P) } , m > 1 

(remember that pm'' = for to < 0). Now one can verify that the unitarity conditions (|70|)-([7l|) are the 
consequencies of the equations ( |87| ) and (|8^). 



6 Photon statistics 

To evaluate the mean number of photons and the statistical properties of the quantum field created in the 
cavity ai t > T we introduce the Hermitian quadrature component operators 

Qm = {am + aL) /\/2, Pm = {am - air) /(*v^)- 
Their variances are defined as 

Um = (g™> - (gm)^ Vm = {Pm) " {Pm)"^ , 

whereas the covariance is given by 

^m "^{PmQm QmPm} {Pm} {Qm} • 

The average values must be calculated in the state defined with respect to the initial operators 6„ (remember 
that we use here the Heisenberg picture). 



6.1 Initial vacuum state 



The vacuum state is defined by means of the relations 6„|0) =0. In this case, C/„i + Vm = 2Mm"''^^ + 1, where 
J\fm°''^'^ is the mean number of photons created from vacuum in the mth mode. It is given by the single sums 



over index n in Eq. (gTl). Initially, C/„(0) = Vm{0) = 1/2, y„j(0) = 0. Using 
wi = 1, we obtain the following expressions for r > 0, 



and assuming for simplicity 



E 



1 



Vm = 



TO \ ^ 1 



(89) 



r„,, = V^Im 



(«)*(«) 

rm r—m 



(90) 



where the coefficients p+^ i, sh ould be taken at the moment T, thus their argument is tt = ^eujiT. Strictly 
speaking, the expressions (p9|)-(po|) have physical meanings at those moments of time T when the wall returns 
to its initial position, i.e. for T = Ntt/[p{1 + S)] with an integer N. Consequently, the argument tt of the 
in (^9|)-(po|) assumes discrete values t'^^ = Ne'K/[2p{l + S)]. One should remember, however, 



coefficients 



that something interesting in our problem happens for the values r ~ 1 (or larger). Then iV ~ e ^ 1, 



£ is so small that tt can be considered as a continuous variable (under 
|l8Sf| ). For this reason, we omit hereafter the subscript T, writing simply 



and the minimal increment At 
the realistic conditions, e < 10^ 
T instead of tt or t^-^\ 

Differentiating the right-hand sides of equations (^) and (Q) with respect to the 'slow time' t, one 
can remove the fraction 1/n with the aid of the recurrence relations ( |87| ) and (|8^). After that, changing if 
necessary the summation index n to n ± p, one can verify that almost all terms in the right-hand sides are 
cancelled, and the infinite series are reduced to the finite sums: 



dUm/dr 
dVm/dT 



p-i 



am 



ip- 



TP. 



in) 
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n=l 

Now one should take into account the structure of the coefficients : they are different from zero provided 
the difference between the upper index n and the lower one m is some multiple of the number p. If m = j +pk 
with J = 1, . . . ,p — 1 and fc = 0, 1, 2, . . ., then only the terms with n = j oi n = p — j survive in the sums 
above. Depending on whether j = p/2 or j ^ p/2, we obtain two different sets of explicit expressions for the 
derivatives of the (co)varianccs. 

1) If m = j + pk but j 7^ p/2 (in particular, for all odd values of p), then 



dUm 

dr 



dVm 
dr 



dr 



0. 



(91) 



In this case Ym = and Um = Vm ~ Mm"'''^ + 1/2. 

2) A different situation happens in the distinguished modes with the numbers fi = p[k + 1/2), k = 
0,1,2,...: 



dU^ 
dr 



-/xRe 



,,(p/2) _ „(P/2) 



dYp/dr = /ilm 



= /iRc 

dr 



.(P/2) 



(P/2) 



(92) 
(93) 



We shall call such modes as "principal" ones; they exist only if p is an even number. In the strict resonance 
case (7 = 0) all the coefficients pjf^^'' are real, so = and dU^/dr < in the whole interval < r < 00, 
resulting in the inequality U^{t) < 1/2, which tells us that the field occurs in the squeezed quantum state. 



6.1.1 Squeezing in the "principal" modes 

(pn-\~n/'2') 

Note that the coefficients depend on the parameter p only through the dependence of the the 

variable k on the product pr: see equations (^9|) or (|l]). Thus to study the squeezing properties of the 
field created due to the NCE, it is sufficient to consider the most important special case of the parametric 
resonance at the double fundamental frequency 2u;i (i.e. p = 2), since the formulae for p > 2 can be obtained 
by a simple rescaling of the 'slow time' (for the 'principal' modes). In this case, only the odd modes can be 
excited from the vacuum, and they do exhibit some squeezing. 

Using equations ( p2[ ) and (|9^) one can immediately find the Taylor expansions of the (co)variances at 
r (assuming (—1)!! = 1): 



V2m+1 



1 

- =F T 

2 ^ 



2m+l 



(2m- 1)!! 



m! 



It 



2m + 1 

(m + ly 



r + 0(r2) 



(94) 



>2m+i = -27(2m+l)r 



2(m+l) 



(2m- 1)!! 



(95) 



We see that the JJ-variances are always less than 1/2 at the initial stage, but the degree of their squeezing 
rapidly decreases with increase of the number m. Note that the dependence on the detuning parameter 7 in 
the short-time limit appears only in terms of the order of r^™^'^ (and higher). 

In the opposite limit r ^ 00 (or k ^ 1), using equations (p^, ( p3| ) and the asymptotics of the Bogoliubov 
coefficients (pq) we obtain constant time derivatives 



dC/2m+i/dr|^_ 
dV2m+i/dr|^_ 

dY2,n+l/dT\^_ 



16a 



7r2(2m + 1) 



16a 



7r2(2m + 1) 
8a 



sm 



cos 



m + 



7r2(2m+ 1) 



sin [(2m + 1) i 



(96) 
(97) 
(98) 
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where (j) = arcsin7. Consequently, all the (co)variances increase with time linearly, giving the constant 
photon generation rate in the 'principal' (odd) modes 



(W2m+l/dT\^_ 



8a 



7r2(2m + !)■ 



(99) 



Equation (99) results in a simple estimation of the mean photon mmrber in the /ith mode at r > 1: M^{t) 
ar I \x. 

Since the covariance Y^^ is different from zero if 7 ^ 0, the initial vacuum state of the field is transformed 



to the correlated quantum state P23, 271 , 272 1. One should remember, however, that the values of C/^, 



and yield the (co)variances of the field quadratures only at the moment t = T (when the wall stopped to 
oscillate). At the subsequent moments of time the quadrature variances exhibit fast oscillations with twice 
the frequency of the mode. For example (omitting the mode index), 

ag{t') = Ucos^icut') + Vsm^{ujt') +Ysm{2ujt'), t' ^t-T. 

Therefore the physical meanings have not the values [/^, V^j and themselves, but rather the minimal 

= Vn values of the quadrature variances during the period of fast oscillations 



u,, and maximal Or, 



73,274| 



} = ^c^M + T v^(c/,,-K<)' + 4y^2) . 



(100) 



Only in the special case of the strict resonane (7 — 0) we have = and Vp_ = V^. In the generic case 
7 ^ 0, all three (co)variances, J7^, V^, and y^, linearly increase with the interaction time T if = t ^ 1, 
due to Eqs. (|9^)-(|9^). Nonetheless, the minimal variance tends to a constant value at t ^ cx). This is 
shown in the subsection 3.2, The examples of explicit time dependences of the coefficients and are 
given in subsection 



6.1.2 Mean photon number 

Differentiating the "vacuum" part of sum ( p7| ) with respect to r and performing the summation over the 
upper index n with the aid of (^7|)-(p8|) (remembering that the coefhcients pm"* arc different from zero 
provided the difference n — m is a multiple of p) one can obtain the formula for the photon generation rate 
from vacuum in each mode (0 < J < P — 1, <Z = 0, 1, 2, . . .) 

d 

l\l \ J _ _V/t( t Tin 1 Ho , ,, 



dr 



2fT(j+TO)Re 
sin(7rj7p)r(g + j/p)T{l + q + ]/p)V{2 - j/p) ^^^,+1 



2pv/r^ ^Y{jlp)T{q+l)T{q + 1) 

xF{q + j/p, ^j/p; 1 + q; K^)F{q + j/p, l-j/p:2 + q; k^) (101) 
We see that there is no photon creation in the modes with numbers p, 2p, .... In the short-time limit, 

r2«+\ r«l. 

In the long-time limit the photon generation rate tends to the constant value (if 7 < 1) 



j+pq 



_d_jr{vac) 

dr 



2ap^ sin^(7rj/p) 
tt'^U +pq) 



1 



0(||ln^ 



apT ^ 1 



(102) 



For g 3> 1 and for a fixed value of k one can simplify the right-hand side of (|101| ) using Stirling's formula 
for the Gamma-functions and the easily verified asymptotical formula 



F{a, b\ c; 2) « (1 — az/c) 



-b 



a, c > 1 . 



In this case 



dr 



MZi - 2pVT^^ sin(7r,7p)r(2-,7p).^ 



2q+l 



7rr(j/p)g2(i-j/p) (1 - k2) 



g > 1. 



(103) 
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In particular, ii S'^{pt)\ ^ 1, then 
d 



(104) 



Comparing ( |102| ) and (104) one can conclude that the number of the effectively excited modes (i.e. the 
modes with a time independent photon generation rate) increases in time exponentially, approximately as 
S^{T)/lnSiT). 

The total number of photons generated from vacuum in all the modes equals 



E 



(105) 



Differentiating ( |105D with respect to r and performing the summation over m with the help of equations 
(p9|)-(|65|) or ( |67| ) one can obtain the formula 



dr 



CO ^ p 

2aRe "^-"^mip- m) p^l'^ (t) pl^'l^r) . 



n 

n— 1 m— 1 



(106) 



Evidently, the right-hand side of this equation equals zero in the "semi- resonance" case p — 1. 

Differentiating equation (106 ) once again over r one can perform the summation over the upper index 
n with the aid of equations (^7)-(|88|) to obtain a closed expression for the second derivative of the total 
number of "vacuum" photons 



dr2 



j^i.ac) ^ 2Re ^ m(p - m) ^t^i'-'Z'^ + vi^-'^H-^^ 



p-i 



2 m{p — m) < mijp — m) 



K ^ / m 



-F 



IP 



; 2; 



in m 



In the short-time limit one obtains 



AA("-) = ip(p2 _ 1), \apT\ « 1 



(107) 



(108) 



In the long-time limit the formulae (|80|), ( |84| ) and X]m=i sin^(7rm/p) — p/2 lead to another simple expression 
(provided p>2) 

_^(™c) ^ 2aV/7r^ a]9T>l, a>0 (109) 

Consequently, the total number of photons created from vacuum due to NSCE increases in time quadratically 
both in the short-time and in the long-time limits (although with different coefficients). 

6.2 Arbitrary initial conditions 

For an arbitrary initial state of the field one can write Um = Um°''^^ + A[/m, where Um°''^^ is given by equation 
( ^) ; similar expressions can be written for Vm and Ym ■ The corrections due to the nonvacuum initial states 
are given by 



A[/„ 



± 



nj 



rm ~ r—rn 



(110) 
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/nj 



rrn r—rn rm I —m 



(6„&j) - {bn){bj 



(111) 



where the average values like 
case of the initial coherent state, 6„|a 



Jj&j) are calculated in the initial state. All the corrections disappear in the 



an I a). If the initial density matrix is diagonal in the Fock basis 



( as h appen s e.g . for the Fock or thermal states) then {b\J)j) = VnSnj {vn > 0), all other average values in 



(110) and (111) being equal to zero. In this case the double sums are reduced to the single ones: 



3^"^ Ay =mV — 



(n) 
P-rn 



(112) 



AF™ 



2m — Im 



(«)*(") 

rm r—m 



(113) 



We see that the initial fluctuations always increase both the variances Um and Vm (for the diagonal density 
matrix). However, asymptotically at r — > oo the corrections are bounded for the physical initial states having 



finite total numbers of photons, because the coefficients pm"* ± 
the summation index n in this limit. For example, if p = 2, then 



and Im 



„(")*„(") 



pt:il\r»i) 



(2«+l) 
P-2' 



'l(T»l) 



2(-l)™ 
7r(2m+ 1) 

2(-l)™ 



(a + z7)™+"+i, 
(a + i7)«-". 



do not depend on 



(114) 



(115) 



7r(2TO + 1) 

and the exponent n disappears in the sums, because \a + i^\ = 1. Thus we have in the 'principal' /i-modes 
(taking p = 2 for the sake of simplicity) 



AC/; 
AVf 
AK 



(oo) 
u 

(oo) 
(oo) 



8Z 



where 



arcsin7, Z — 



2sin^ {H(j)/2) 
2cos2 {fi(j)/2) 
— sin {ficf)) 



(116) 



^ 2fc + 1' 



The expressions in (116) are very similar to those in equations (|96|)-(|98|). The consequence of Eq. (116) is 
the important result that in the limit r — > oo the minimal variance does not depend on the initial state 
of the field inside the cavity, provided the initial density matrix was diagonal in the Fock basis. Indeed, in 
this case, combining the equations (|96|)-(|9^) and (|ll6| ), we can write the variances at t ^ 1 as (we omit the 
subscript p) 

U{t)\ / 2Fsin2(x/2) + / 
Vir) = 2Fcos2(x/2)+g I , F 
Y{t) I \ -Fainx + h 



[OT + . 



X 



M<^. (117) 

The corrections /, g and h can be found by integrating equations (|92|)-(|93|), therefore they do not depend 
on the initial state. At t — )■ oo these corrections tend to finite limits, so they are much smaller than F. 
Evidently, U + V = 2F + f + g, whereas 

{U - Vf + = 4^2 ^ 4^[(^^ _ J) cosx - 2/isinx] + (./ - gf + 4/1^. 
For F ^ f,g,h we have 



^/{U -V)^ +4Y^ = 2F + {g ~ f)cosx-2hsinx + 0{l/F) 
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so the minimal variance u{t) ( |100| ) tends to the finite Umit 

u(oo) = /cos2(x/2) +5sin2(x/2) + hsmx, 



(118) 



which does not depend on Z, i.e. on the initial state. 

The correction to the mean number of photons in the m-th mode for the "diagonal" initial distributions 
is given by the sum 



n 



+ 



An) 



It tends to the limit AJ\f^°°^ = 8Z / (tt^^). 

The total number of photons in all the modes equals Af = A/"'""^-* + A/"*^^"^' , where 



ik 



(119) 



(to obtain this formula one should use the identity (|2J)). Differentiating (119) with respect to r and 
performing the summation over m with the help of equations (|2|), (|3^), (g), ^), or (|7|), one can obtain 
the formula 



dT 



= 2cr > > m{p - m) 

•/nk 



7i,k—l 



2crRe ^ "— m(p — m) 



(n)* {k) (k) (n)* 

P—m Pp—m ' P—mPp—m 



P 



.,k=l 



' nk 



m—l 

An) 



' in) (k) ^ (n) Ak) ■ 
r —ml rn—p ' rm I p — m 



(120) 



Using equation ( |120|) and replacing the coefficients pm by their asymptotical values (^5|) one can obtain the 
expression 



= — ^ ^ sm (^to/p) ^ 



,n+fe 



— 1 



{(^L+p«Wp/=>(a + «7) 



k — n 



(jRe 



^ ^(m + pn) (m + pk) 
{bm-}-pn^m-\- pk){a + i-f) 



k+n+l 



(121) 



which holds provided apr 3> 1 and a > 0. For the physical initial states the sum in the right-hand side of 
(121) is finite. This is obvious if a finite number of modes was excited initially. But even if the cavity was 
initially in a high-temperature thermal state, so that {h\fik) — ^nk'd /n, (bnbk) = 0, the sum over n, k yields a 
finite value Q X^J^o (™ +pn)~^. Consequently, the total number of "nonvacuum" photons increases in time 
linearly at apr ^ 1, whereas the total number of quanta generated from vacuum increases quadratically in 
the long time limit. 



6.3 The "principal resonance" {p = 2) 

Many formulas obtained above can be simplified in the special case p = 2. In this case there are two subsets 
of nonzero Bogoliubov coefficients. The first one consists of the coefficients with even upper and lower indices 



^2k'^ which are reduced to the coefficients S}^^ of the "semi-resonance" case. However, since r]^^'^ = 0, this 
subset does not contribute to the generation of new photons. The second subset is formed by the "odd" 
coefficients which can be written as [n = k(2t)] 



f(2n+l) 
S2m+1 



r(n + 3/2)K"-'"A™+"+i 
r(m + 3/2)r(l + n-TO) 

(n + 1/2 , — m — 1/2; 1 + n — m; k^) , n > m 



(122) 
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(2„+i) ^ (-l)'"-"r (m + 1/2) k"'-"A"'+"+i 
r(n+l/2)r(l + m-n) 
X (to + 1/2 , -n - 1/2 ; 1 + m - n; K^) , m > 



(123) 



(2„+l) _ (-l)^-lr (fc + 1/2) r (n + 3/2) ^n+k+l^n-k 

'^2'=+! " (2 + n + fc) 

X F(7i + l/2, fc + 1/2; 2 + n + fc; K^) . 



(124) 



It is known |275| that the hypergeometric function F{a, b; c; z) with 'half-integral' parameters a, 6 and an 
integral parameter c can be expressed in terms of the complete elliptic integrals 



K(«) = 



7r/2 



da 



TT /I 1 



\/l — sin^ a 



TT / 1 1 \ 

E(k) = y daVl - sin^a = f (^^ 2 ' 2 ' ' j " 



2^ 2' 2^^^'^^ 



In particular, 



^ ^a(k)E(k), r;!'^ = — rK2K(K) _ E(ac) 

TT TTK 



P3 



P-3 = 



^'' = ^[(1-2^^^)E(«.)"-'K(.)] 



37rK 
2 



[(2-k2)e(k)-2k2K(k)] 



where 



and A(k) was defined in (pO|). 



KEE - = [l + S'2(2r; 



1-1/2 



The general structure of the coefficients p^^^ (we confine ourselves to the case p = 2) is as follows 

P^lUl = ^^^^ [fm («^) E(«) + khra («^) K(..)] 



(125) 

(126) 

(127) 
(128) 
(129) 

(130) 
(131) 



[rra (k') E(k) + K^s,^ (^2) K(k)] 



(132) 



7rK'"+lA™(K) 

where fmjx), gm{x), r„i{x), s„i{x) are the polynomials of the degree to which can be found from the recurrence 
relations (|67|). 

The photon generation rate from vacuum in the principal cavity mode (m ~ 1) reads 



(vac) 



dr 



-2Re 



E(k) [E(k) - k2K(k)] 



(133) 



The average number of photons in the first mode can be obtained by integrating this equation. It is convenient 
to integrate with respect to the variable k, taking into account the relation (for p — 2, for instance) 



dK = 2(ik^AT, /3 = ReA = 

For example, in the case of the quadrature variance t/i we arrive at the equation 

AUi 2 



dK 



,2^2,2^ (1 - 27 V) + 1 - 2/?k] E2(.) 
-2^2(1 - /3k)E(k)K(k) + K^K^in)] . 



(134) 
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Let us consider first the case 7 = 0, when /3 = 1. Taking into account the differentiation rules ||27i 

dK{K) _ E(k) K(k) dE(K) _ E(k) - K(k) 



dK 



dK 



(135) 



we may suppose that the factor in the denominator of the right-hand side of equation (134) comes from 
the derivative dK/d/t. Thus it is natural to look for the solution in the form 



Ui^ — [A{k)K^{k) + B(k)K(k)E(k) + C(k)E^(k)] 



(136) 



where A{k), B{k) and C(k) are some polynomials of k. Putting the expression (136) into equation (134) 
we obtain a set of coupled equations for the unknown functions A,B,C. Writing A{k) =00 + ^i(k), 
B{k) = 60 + Bi[k), C{k) = Co + Ci(k) we determine the constant coefficients ao, 60 and cq by putting 
K = in that equations. Then we obtain new equations for the functions Ai(k), Bi{k) and Ci{k) and 
repeat the procedure. After a few steps we arrive at the equations which have obvious trivial solutions 
An = Bn = C„ = 0. This confirms our hypothesis on the polynomial structure of the functions A{k), B{k) 
and C(k) and gives the final answer. The equations for the variances J7^, V^, etc. with > 3 can be integrated 
in the same manner, the only difference is that one should write k'^ instead of k in the denominator of the 
expression like (136). In the generic case 7 ^ we notice that the factor P can appear in the denominator 
of the expression (134) as a result of differentiating the function (3{k), since d/3/dK — — 7^k//3. Therefore we 
split each function, A,B,C in the '/9-even' and '/3-odd' parts like A = Ae{K) + /3{k)Ao{k). The equations 
for the 'even' and 'odd' coefficients turn out independent, and we solve them using the procedure described 
above. 

The results of the integrations are as follows, 



AA(-)(«) = ^Kin) [2E(«:) - ii'K{n)] - 1. 

^ [k^{(3 - k)K^{k) - 2{f3 - k)K(k)E(k) + PE'^{k)] , 
[2(/3 + k)K(k)E(k) - k^ip + k)K^{k) ~ fJE^in)] , 



(137) 



2 
2 



Yi^^ \k^K^(K) - 2K(k)E(k) + E^U)] 



Making the transformation [269, 27C] 

'1 - V 



K 



1 



-K(k), E 



1 - ft 

1 + K 



E(k) + kK{K) 



1 



one can rewrite formulae (127) and (137) in the form found for the first tim e in [189| (in the special case of 
7 = 0). Using the asymptotical expansions of the elliptic integrals at k — > 1 |276] 



one can obtain the formula 

^(vac) 



E(k) « 
(r » 1) = -^T - 



In - + - I In - - 1 



4 1 
K 4 
1 
2 



In- 



- 2 

K 



- + {re-'^"^) , a > 



(138) 



In the special case of 7 = 1 one can obtain the expansion 



{vac) 



r > 1) 



■ In T ■ 



12 



■In 2 



+ 0(r-2) 
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If 7 > 1, the number of photons in the principal mode oscillates with the period tt/ (25). For 7 3> 1 one can 
write K « sin(2ar)/a, i.e. \k\ <^ 1. In this case 



j^ivac) _ ^ p. sin^(2aT) ^ ^ 
^ 4 ~ 

Equation ( |lOO|) yields the minimal and maximal invariant variances 

2 

wi = — \k^(l - k)K'^(k) - 2(1 - k)K(k)E(k) + E2(k)1 , (139) 

t;i = ^ [2(1 + k)K(k)E(k) - k2(1 + k)K2(k) - E2(k)1 , (140) 

which depend on the detuning parameter 7 only implicitly, through the dependence on 7 of the function 
k(t). In the short time limit t ^ 1 (then k w 2r) wc obtain, using the Taylor expansions of the complete 



elliptic integrals, mi = ^ — t + + • • • and wi = ^ + t + + • • • in accordance with [124|. More precisely, 



The minimal variance ui monotonously decreases from the value 1/2 at i = to the constant asymptotical 
value 2/7r^ at r ^ 1, confirming qualitatively the evaluations of |107, 11C| and giving almost 50% squeezing 



in the initial vacuum state. The variance of the conjugate quadrature monotonously increases, and for 
T ^ 1 it becomes practically linear function of time: ui(r 3> 1) ~ 16r/7r^. The asymptotical minimal value 
ui{t — 00) does not depend on 7 provided 7 < 1 (only the rate of reaching this asymptotical value decreases 
with 7 as \/l — 7^). In the strongly detuned case, 7 > 1, the minimal variance oscillates as a function of t 
(being always greater than 2/7r^), since in this case the function k{t) oscillates between and 7~^. 

The minimal variance does not go to zero when r — > 00 due to the strong intermode interaction, which 
results in a high degree of quantum mixing for each mode. Since the state originating from the initial vacuum 
state belongs to the class of Gaussian states (see the next subsection), the quantum 'purity' Xm = Trp^ 
of the mth field mode (described by means of the density matrix p^) can be expressed in terms of the 

(co)variances as [ ^77| Xm = [4 (CAnK™ — ^m)] Using equations (117)-(118) one can check that for 



T 3> 1, UV — = 2Fu{oo) + 0(1) ~ T. Consequently, the purity factor x asymptotically goes to zero as 
^-i/2_ ppj, instance, for m = 1 we have (writing simply K and E instead of K(k) and E(k)) 

XI = [4KE3 + m^K'E - ea^K^E^ - E^ - ^^K^] (141) 

The initial dependence on k is rather weak: x('* ^ 1) = 1 — i^'^^ + ■ ■ ■• But when k ^ 1, x rapidly goes to 
zero: x(k < 1) w (S/vr^) [In (4/k)]"^^^ with dx/dn -> ~oo. 

The expressions for the variances in the modes with numbers /i > 3 are rather involved. Here we give 
only one explicit example - the variance C/3 for 7 = 0: 

+ (1 - k) (4k3 - Uk^ - 20k ~ 8) K(/c)E(k) 

+ (4k'^ + 6k^-k2 + 6k + 4)e2(k)] (142) 



The Taylor expansion of the right-hand side of ( |142| ) coincides with the expansion (|9J). The asymptotical 
value at r — !■ cxD equals U3{k = 1) = 38/(97r^) « 0.43. We see that the squeezing rapidly disappears 



with increase of the mode number /i. The variance V3 can be obtained from (142) by means of a simple 



substitution k — > — k. Therefore the mean number of photons in the third mode is given by 



■^3 = [(3^' - 2) K (2E - ^^K) + 2 (l + k') E^] - - . (143) 
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The second derivative of the total mean number of photons created from vacuum in aU the modes (107) 
takes the form 

- ^ \^'^' - 2«^KE + (1 + - 27^«^) E^] . (144) 

Integrating this equ ation with the account of the condition AM''^"''^'' /dT = at r = [which is a trivial 
consequence of Eq. (106)], we obtain very simple expression 

In the limiting cases this formula yields 

J^(vac)(^^ < 1) W r2 

AA('"^=)(r> 1) = 8aV/7r2 + 0(r), a > 0. 

If 7 ^ 1, then |k| <C 1, but 7^k^ « sin^(2aT) ~ 0(1). In this case the Taylor expansion of the second 
derivative yields A/''^*"^'^' = 2 cos(4aT)+0(7~^). Integrating this equation with account of the initial conditions 

jif{vac)(^Q^^ = 7V(^"'=)(0) = one obtains M^'"'"'^ « M['''"'^ « sin^ {2hT) / {Ah"^) . 

6.4 Photon distribution 

Now let us turn to the photon distribution function (PDF) /(n) = (n|pm(^)|n), where |n) is the multimode 
Fock state, n = (ni,n2, . . .)j ^-nd pm{t) is the time-dependent density matrix of the mth field mode in the 
Schrodinger picture. Note that all the calculations in the preceding sections were performed in the framework 
of the Heisenberg picture. Nonetheless, the available information is sufhcient to calculate the PDF for the 
special (but very important) class of Gaussian initial states (defined as the states whose density matrices, 
or wave functions, or Wigner functions, are described by some Gaussian exponentials). This class includes 
coherent, squeezed and thermal states; in particular, it includes the vacuum state. 



The solution is based on two key points. The first one is the statement [170 



176 



182 



185] that the 



field evolution in a cavity with moving boundaries can be described not only in the Heisenberg picture, 
but, equivalently, in the framework of the Schrodinger picture, with a quadratic multidimensional time- 
dependent Hamiltonian. The second key point is the fact | 278| , 279 1 that the evolution governed by quadratic 
Hamiltonians transforms any Gaussian state to another Gaussian state. Knowing these facts, it remains to 
take into account that the photon distribution function of any Gaussian state is determined completely by 



the average values of quadratures and by their variances [280, 281], which obviously do not depend on the 
quantum mechanical representation. 

In the most compact form the information on the photon distribution f{n) in some mode (we suppress 
here the mode index) is contained in the generating function 



G{z) ^ J2 /(")^" 



For a generic one-mode Gaussian state it can be expressed as [192, 280, 281| 



where 



G(z) = [c;(z)]-V2exp(- 



zgi - z^92 



90 



g{z) = ^ [(1 + zf + 4{uv- y2) (1 _ 2)2 + 2{u + v){i- z^)] 

D = 1 + 2{U + V) +4 {UV - Y^) = 4^(0) 
go = {pf{2U + 1) + {qf{2V + 1) - 4{p){q)Y 



(146) 
(147) 



9i = 



2{p)^ iU^+Y^ 



U 



2{qfiV^ + Y^ 



V- 



~4{p){q)Y{U + V + l) 
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92 - 2(p)2 {u^ + y2 _ 1 j ^ 2(g)2 {v^ + y2 _ 1 ^ _ 4^^^ i^^^y{U + V). 
In the generic case f{n) is related to the two-dimensional 'diagonal' Hermite polynomials 

f[n)^^Hi^Hx,x*) (148) 

where 

To = f{Q) = 2D-^/^e^p[-g^lD) 

_ V2{{2V ~ l){q) - 2Y{p) 2U){p) + 2Y{q)]} 

^ 2{U + V)-4:{UV-Y'^)-1 

and 2x2 symmetric matrix TZ has the elements 

7^ll = 7^^2 = ^{v -u- 2iY), 7^l2 - 7^2l - [i - 4 (c/y - f^)] 



The two-dimensional Hermite polynomials are defined via the expansion |269[ 

exp --a7^a + a7^xU ^ (149) 



m.n— 



where x = (si, X2), a = (ai, 02). The properties of these polynomials were studied in |281, ^82 |. In particular, 
they can be expressed as finite sums of the products of th e usual (one-dimensional) Hermite polynomials. 
The corresponding formula for the probabilities reads |280| 



where 



A = yJ{U -Vf +AY'^, S = A{UV -Y"^) -I 
^ ^ {2V + l)(g) - 2Y{p) + 2U){p) - 2Y{q)] 

[2D{V ~U -2iY)f''^ 

If {p) = (q) = 0, then the generating function ( |146D is reduced to [t/(z)]^^/^, i.e. it has the same structure 
as the known generating function of the Legendre polynomials P„ (x) . In this case, we have the following 
expression for the photon distribution in the mth field mode: 

r ( , 2[(2u,„-l)(2z;^-l)]"/^ ^ / Au^v^ - 1 ^ 
[(2u„ + l)(2z;„ + l)]("+i)/2 "1^V(4<-1)(4«,1-1); 

It depends only on the invariant minimal and maximal variances Um and Vm- Note that the argument of the 
polynomial in (151) is always outside the 'traditional' interval (—1,1) (in particular, this argument is pure 
imaginary if 2um < 1), being exactly equal to 1 for the 'non-principal' modes with Um = Vm = ■N'm + 5, 
when formula (151) transforms to the time-dependent Planck's distribution 

^ ' [AA™(r) + iri 

Only for the 'principal' /i-modes the spectrum of photons is different from Planck's one due to the squeezing 
effect. 

The function (151) can be simplified in the long-time limit t ^ 1, when the average number of created 
photons M = n K, {V + U)/2 exceeds 1. Then the mean-square fluctuation of the photon number has the 
same order of magnitude as the mean photon number itself, y^ovT ~ V^A/", and the most significant part 
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of the spectrum corre spon ds to the values n ^ 1. Using the Laplace-Heine asymptotical formula for the 
Legendre polynomial | 283 | 

V2^(^2_ 1)1/4 



one can simplify (151) for the fixed values of the invariant variances u and v as 



fin) 



y^Trn{v 



2v-l 
2v + l 



n+l/2 



(152) 



provided the positive difference v — u is not too small. Another approximate formula can be used if u 2> 1 
but u 1: 



The first and second derivatives of the generating function (146) at z = 1 yield the first two moments of 
the photon distribution (hereafter we suppress subscript to) 



n = i(u + w - 1), CT„ = n2 - {nf = i {2u'^ + 2v'^ - l) 



which result in the Mandel parameter | 284 | 



= (Tn/n- 1 



+ v'^ -u-v + l/2 

U + V — 1 



(154) 



(155) 



This parameter appears positive for all values of r, so the photon statistics is super-Poissonian, with strong 
bunching of photons (the pair creation of photons in the NSCE was discussed in |4|,|l98[|o|,|213)). In 
particular, 

Q2™+i(t ^ 0) « [(to + 1)(2to - 1)!!/to!]' t2'"/(2to + 1), Qi(0) = 1, 
whereas Qm ~ Kn(T) 1 for t ^ 1 (if 7 ^ 1). 

7 Energy and formation of packets 
7.1 Energy density 

The mean value of the energy density operator in one space dimension 



W{x,t) 



(dA/dt^'^ + (dA/dxy 



(156) 



at i > T equals (hereafter we assume — 1, i.e. oji ~ tt) 



W{x,t)=Ti ^ cos[7r(TO + j>]Re (a^a^^-^f^+J)*' 

m.j — 1 \, 

+ icos[^(TO- j>] [(aj„a,)e''"(™-^'*' + (a„aJ)e-"(™-^)*' 
2 L -' 



(157) 



where the quantum mechanical averaging (• • •) is performed over the initial state of the field (the Heisenberg 
picture) and t' = t + ST. 
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7.1.1 Regularization and Casimir's energy 



Due to the commutation relations 
part with m = j: 

oo 

yirivac) ^ (^^2) mexp{-iTTmt' + iirmt') 



Smj the series (157) contains the vacuum divergent 'diagonal' 



(158) 



The recipe how to regularize this divergence was given in ||81| . One should write the first term in the 
argument of the exponential in (158) as it stands, but to replace t' in the second term by t' + iri, r] > (this 
is the so-called 'point-splitting method'). Then the sum becomes convergent, giving 

oo 

W^(™'=)(r;) = (7r/2) ^ me"™'^'' = (tt/S) [sinh(7rry/2)]"^ . 

m— 1 

The Taylor expansion of this function reads W^'''"''> {-q) = i2nr]'^)-'^ - 7r/24 -f- ©(jy^). According to one 
should remove the divergent term (27r77^)~^ and after that proceed to the limit rj ^ 0. This limit value gives 
us the known expression for the one-dimensional negative vacuum Casimir energy ^,|rTW20| (which does 
not depend on the coordinate x in the case involved) 



yiricas) ^ _7r/24 (or 



nhc 
2412 



in the dimensional units) 



(159) 



Extracting this vacuum energy from W we arrive at the expression 



W = =7r ^ y^Re (a„aj)cos[7r(m-Hj)a;]e-'^(™+^')*' 

m J — 1 

+ {alaj) cos[7r(m - j>]e"(™~J)*' 



(160) 



The same expression ( |160| ) can be obtained if one calculates the mean value of the normally ordered (with 
respect to the operators aj^ and a„) counterpart of the operator (156) (cf. [ Il64| ). Then the total energy 
(without the vacuum part) assumes the usual form 



£ = 



W{x, t)dx = ^ w„ {alfin) 



(161) 



which justifies the choice of the normalization in ( |12D and (|2l|). 

Since the initial quantum state was defined with respect to the 'in'-operators 6j, and 6„, we must express 
the 'out' operators and am in terms of 6^ and 6„ by means of formula (|2^). Thus we arrive at the 
expression containing a combination of the mean values (6„6fc), (6]j6^), {bU>k) and (bnb\}i calculated in the 
initial quantum state. For the initial vacuum state defined according to the relations fe„|0) — 0, n = 1, 2, . . ., 
the only nonzero mean values are (J)ni>\) — Snk- Then (|16C|) is transformed into the triple sum 



Woix,t) 



{ cos[^(m - j)x]e"("-^')*'pLlpL'f 



^ n 

n.m,j — l 

- cos[7r(m + j>]e"('"+^)*V^lpf 



(162) 



Evidently, Wo{x,t) = for t < 0. For an arbitrary initial state the energy density can be written as a sum 
of the 'vacuum' and 'nonvacuum' contributions 



1 



=Re 



(163) 



n,/c— 1 
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where 



^(nk) ^ ^ cos[7r(m + j)a 

mj"=l L 



^i7r(m+j)t' (n) (fc) 



iw{m+j)t' (n) (k) 



■ COS 



(164) 



j^{nk) ^ ^ ^^-1 cos[7r(m - 

mj — 1 



I i7r(m-j)t' „(")* ^C^') 



„-iTT(m~j)t' (n)* (k) 



■ cos[7r(m + j)x\ 



^iTT(m+j)t' M* (k) i7rim+j)t' (71)* (k) 



(165) 



Making the change of the summation index j — > — j in the first term of (|162| ) we can write 



OO OO CO 



W^{x, t) = -TrRe ^ E E V ''^'t^^™ + J>]e"(™+-'")*'pLlp5"^*. 

71—1 m—1 j— — oo 



Similarly, changing the indices m —m or j — > —j in (164) and ( |165| ) we can reduce 4 sums with apparently 
different sumniands and the indices running from 1 to oo to the unified sums whose two indices run from 

—OO to oo: 



j^{nk) _ ^ mj cos[7r(m + j)a;]e 



^("i+j)t' M (k) 
rm rj 



mj — — OO 
oc 



^(nk) ^ ^ mj cos[7r(m — j)x]e 



(m-j)t' Jn)* (k) 
rm rj 



Now, replacing the cosine function by the sum of two imaginary exponentials we see that W{x, t) is actually 
the sum of two identical functions of the light-cone variables: 



W{x,t)^ -[F{u]t) + F{v:t)], u = t' + x, v^t'-x, 



where 



F = F^+ E 



' nk 



=Re 



{bJk}F^''''^ + {blh}F^"''^ 



n,k—l 

OO OO OO 



Fo{u;r)^-ReY^Yl E ^e-('"+^>pLl(r)p;")*(r) 



n— 1 m—1 j^ — OO 



(166) 

(167) 
(168) 



OO OC 



F^-'\u;r) = E E "^Je""<"+^'^>l:HT)pf (r) (169) 

ni— — OO J — — OO 

OO OC 

F("'=)(ii;T) = ^J^""^"'''''^"Pn?*i^)pf\^) (170) 



m— — OO J — — 00 



The extra argument r in the above expressions is introduced in order to emphasize that the energy density 
depends not only on the value of the current time variable t (which must satisfy the condition t > (1 + d)T), 
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but also on the moment of time T when the wall stopped to move. It is worth mentioning that the variables 
t and T are independent, as well as u and r or i; and r. 

Evidently, the double sums (16£) and ( |l70| ) are factorizcd to the products of independent sums over ni 
and j: 



m— — oo 



-i-rrmu (n) ( \ _ ^ ^ 



E 



(171) 



The last sum in (171) can be easily expressed in terms of the generating function ( p'4[ ) if one writes n — j + kp, 
m — j + Ip and z — exp(— iTrpu). Thus we obtain 



G(")(w;t) = 



nz [zg(jiT) + aS{pT)] 
[g*{pT) + zaS{pT)] 



where 



f{u;K) ^ \g*{pT) + zaS(pT)\ = 



2— exp( — 27rpu) 



1 -K^ 



[1 + + 2(TK COs(p7r7i — if)] 



2 ' 



^ ^ Zff(pT) + cr5(pr) ^ gj(2>p-7rp«) 1 + crKexp[i(7rpu - (/?)] 



g*{pT) + z<tS{pt) 

Sjpr) 
V^ + SHptY 



1 + (7Kexp[i{ip — T:pu)] ' 



(172) 

(173) 
(174) 



exp 



(iif) = yl — 7^«? + . 



In the 'vacuum' contribution (16S) we have some asymmetry between the indices m and j, since m runs 
from 1 to oo, whe reas j runs from — oo to oo. This asymmetry can be eliminated if one differentiates both 
sides of equation ( 168 ) with respect to the independent variable r at a fixed value of u and performs the 
summation over the superscript n with the aid of the recurrence relations (^^ and ( |88| ) . It is easy to verify 
that all the summands with n > p are canceled, so the infinite series over n can be reduced to the finite sum 
from 1 to (p — 1): 



dr 



p— 1 OO OO 

= -aRe^ ^ J2 rnje 

71—1 7n — l j — — oo 



i7T{m-\-j)u 



[pLl(r)pL7")(r)+p(r")*(r)p(")*(r) 



(175) 



Making the change of summation indices m 



-j, n ^ p—n in the first product inside the square 



brackets one can reduce two sums in the right-hand side of (175) to the single series where both the indices 
m and j run from — oo to oo. Moreover, the sums over m and j become completely independent, giving rise 
to the equation 

^^^^ = -aRe J2 G^"^ [u; t)G(^'-") {u; r), (176) 

n— 1 

where G'^^^\u]t) is given by (171). Due to ( |172| ) the sum in the right-hand side of ( |176| ) is reduced to the 
sum 'YTn=i ^(P ^ ^) = ^p{p'^ ^ !)• Introducing the variable rj — exp(2apT) we obtain the explicit expression 



dFo{u;T]) 
drj 



(p2 - l)a^j^ [jfil + Q,-f/3)+a-/3-l] 
12 [?72(1 + « + /?)- 277(72 + /?) + 1 + /? - a]' 



(177) 



where a — cracos(p7rM) and (3 ~ (77 sin(p7ru). Integrating (177) with the initial condition i^o = at r = 
(or ?7 = 1) we arrive after some algebra at the simple expression 



Fo{u- k)^B [f{u; B={p^~ l) /24 , 



(178) 
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where function /(u; k) is given by (173). Finally we obtain the following expression for the function F{u] t) 
defined by equation ( |l66[ ): 



F=-B + fiu; t)\b+ J2 ^I^e \{bJk)A'^ + {bik)A 



,fc=i 



(179) 



In the special case of the initial states whose density matrix is diagonal in the Fock basis, so that {blfik) = 
t^n^nk and (bnbk) = (for example, the Fock or thermal states; Vn is the mean number of quanta in the nth 
mode), the sum in ( |179| ) is proportional to the initial total energy £q in all the modes (above the Casimir 
level): 

Fid^ag) ^) = _s + f (u; t) [B + A^o] , (180) 



n=l 



7.2 Packet formation 

Now let us analyse the expressions for the energy density obtained above. For the initial vacuum state we 
see immediately from equations (173) and (178) that in the generic case the function Wo{x,t) with the fixed 
value of the 'fast time' t has p peaks in the interval < a; < 1, whose positions are determined by the 
equations crcos(p7ru — ip) — —1 and a cos{p'Kv — ip) — —1. Obviously, for t > T the energy density is a 
periodic function of the time variable with the period = 1 if p is an even number and A< = 2 if p is 
odd. 

For the even values of the resonance multiplicity p we have p/2 peaks moving (with the light speed) in 
the positive direction and p/2 peaks moving in the negative direction. If p is odd, then the numbers of peaks 
of each kind differ by 1 . All the peaks have the same height 



lyi— ) = 27rS«/(« - If = ^B {e'P- - l) 



(181) 



(in this section the expressions containing r are related to the special case of the strict resonance 7 = 0), 
excepting some distinguished instants of time when two peaks moving in the opposite directions merge, 
forming a peak with double the height. 

If K — > 1 (i.e. r > 1 and 7 < 1), then the energy density can be approximated in the vicinity of each 
peak by the Lorentz-like distribution 



W, 



(vac) 



1 + {2Sx/A 



2 ' 



(182) 



where the width 



Ai 



2 1 



/4 



1/4) 



-2pT 



PTT y/K pTT 

of each peak is defined as the double distance between the position of the maximum and the point where 
the energy density decreases 4 times. One can introduce also the 'energy width' of each peak by means of 



the relation W„ 



£{t)/p. For K - 



1 we obtain 



(1 - n)/{2^p) « (ttp) 



V2{ttp)- 



nearby the peaks the 'dynam- 



Except for narrow regions of the length A+ k, {ttp)^^^/\ — 
ical' energy density is less than its initial vacuum value, in agreement with the results of |112, 164, 165, 168| 
obtained in the framework of different approaches. The minimum values of Wq far off the peaks are given 
by (taking into account the contributions of both the functions Fq{u) and Fq{v)) 



W„ 



-AttB 



(183) 
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If K ^ 1, Wmin —T^ {p^ — l) /24. Ad ding to this expression the initial Casimir energy ( |159| ) we obtain the 



total asymptotical minimum value (cf. |164]) 



W^it = - V/24. (184) 

For an arbitrary initial state the energy density has, besides the 'vacuum' part, the additional terms given 
in equation ( |17S| ). Since these terms are proportional to the same functions /(u; k) or f{v; k) which determine 
the structure of the 'vacuum' part, the positions of the peaks are not changed (remember that |A| = 1). For 
the initial states with diagonal density matrices in the Fock basis (in particular, for the thermal states) all 
the peaks still have equal heights, increased by the quantity AWmax — ^^oi^ + '^)^/(l — '«)^: compared with 
the vacuum case. However, the asymptotical minimal value of the energy density at k ^ 1 does not depend 



on the initial state, being given by formula (184) in all the cases. 

If the initial density matrix in the Fock basis has nonzero off-diagonal elements (as happens, in particular, 
for any pure state different from the Fock one, e.g., for the coherent states), different terms in the sum ( |179| ) 
can interfere. As a consequence, the peaks acquire some kind of 'fine structure'. For example, if only the 
first mode was excited initially in the coherent state |a), a = |a|exp(i'0), then for p = 2 and 7 = (the 
strict resonance) we have 



, N „ (l — K^)^ sin(z + -0) + sin(z — ?/;)l'^ 



[(l-K)2+4KSin2 z] 

where z = tt (u — it*) and is the position of the 'vacuum' peak determined above, liip = tt/2, then we have 
the high maximum AVF^y2^ = TT\a\'^ {1 + k)'^ / {1 — k)'^ at z = 0. However, iitp = 0, then instead of a maximum 
we have the minimum AW = at the same point z = 0, and the peak is split in two symmetric humps with 
equal maximal heights AWq™"'^ = ([1 + k\'^/27k) AW^ji^ located at the points sinz = ±(1 - k)/V^. In 
the intermediate case < i/' < 7r/2 asymmetric forms of the peaks are observed. If p > 2 or several modes 



were excited initially, the interference between different terms in (179) can result in different heights of the 
peaks and more complicated 'fine structures' (provided (6„fefc) ^ for some n and k). 

If the detuning 7 is different from zero, then some deformations of the form of peaks are observed, 
although the maximal heights are still of the same order of magnitude as in the case 7 = 0, as far as 7 < 1. 
But if the detuning exceeds the critical value 7 = 1, the energy becomes an oscillating function of the 'slow 
time' r, the amplitude of oscillations being proportional approximately to (7^ — l) [191,092]. The peaks 
become rather wide and low, since the parameter n is limited by the inequality k < 7~Mi 7 > 1. The 



illustrations can be found in |193] 



Since the components of the energy- momentum tensor Tqo and Tn are given by similar expressions in 
the case of single space dimension, the force acting on each wall has the same time dependence as the energy 
density at the points a; = and x = I. For the most part of time during the period of field oscillations 2Lo / c 
(where Lq is the distance between the walls at rest) this force is negative, being less than the static Casimir 
force, with the maximal amplification coefficient p^. However, the average value of the force over the period 
is positive due to the creation of real photons inside the cavity. 

If the walls possess some small transmission coefficient, then a small part of radiation accumulated inside 
the cavity can leave it. In this case one could observe sharp pulses of radiation outside the cavity |210[ , 
whose amplitudes must be proportional to the heights of the peaks inside the cavity multiplied by the 
small transmission coefficient. The intensity of these pulses can be significantly increased, if the initial state 



is different from vacuum and possesses sufficient energy, like thermal states 210 , 211| or coherent states. 
However, to describe the form of the pulses exactly it is necessary to develop a more general theory which 
would take into account the boundary conditions corresponding to the partially transmitting walls (because 
the nonzero transmission coefficient can change significantly the pulse shape, just as the nonzero detuning 
deformed the form of packets in the examples considered above). 

7.3 Total energy 



The total energy (161) of the field inside the cavity (above the initial Casimir level) can be o btain ed by 



integrating the density W{x) (|166| ) over x. The contribution of the vacuum (fimction Fq in ( 167 )) and 
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diagonal' terms (given by the partial sum in ( |167D over n — k) can be calculated with the aid of the formula 



dx 



{a + bcosxY (a2_t2)i/2- 

To find the contribution of 'nondiagonal' terms (n ^ k) it is convenient to replace the integration over x by 
the integration in the complex z-plane (z = exp[— ?7rpu] or z = exp[— zTrpu]) over the circle \z\ = 1. One can 
check that this circle is passed p times when x goes from to 1 (if one takes into account both the 'm'- and 
'ii'-contributions). It turns out that the integrals of the 'non-diagonal' terms are different from zero if and 
only if the corresponding integrands in the z-plane have simple poles inside the circle \z\ = 1. This happens 
only when k + n = p in the first term inside the square brackets in (179) and k — n = p in the second term 
inside the same brackets. 

Finally we obtain a simple expression 



Sir) =£o + 2S\pt) 



So 



where 



-nB+^lmig)\-^S{2pT)Rcig), 



p-i 



g^2T:J2 Vn{n+p){bl 



(185) 



(186) 



n=l 



{if p — 1, the last sum in ( |l86[) should be replaced by zero). Formula (185) was found for the first time in a 
different way in | 191 |. In particular, for the initial vacuum state of field we have 



1 



12a2 



sinh^ {par) 



(187) 



The total energy increases exponentially at r — > oo, provided 7 < 1. In the special case 7 = such 
asymptotical behaviour of the total energy was obtained also in the frameworks of other approaches in 



|164-166, 177|. Here we have found the explicit dependence of the total energy on time in the whole interval 
< T < 00, as well as a nontrivial dependence on the initial state of field, which is contained in the constant 
parameter g. This parameter is equal to zero for initial Fock or thermal states of the field. However, in a 
generic case g 0, and it can affect significantly the total energy, if £{0) 3> 1. Consider, for example, the 
case p = 2. If initially the first mode (n = 1) was in the coherent state |a) with a — |a|e*"^, \a\ ^ 1, and all 



other modes were not excited, then £(0) 



% so for r ^ 1 and 7 = (exact resonance) we have 



£{t 3> 1) « j\a\'^e [2 — cos(2^)]. The maximal value of the energy in this case is three times bigger than 
the minimal one, depending on the phase (j). 

According to (185), the initial stage of the evolution does not depend on the detuning parameter 7 for 
all states which yield Im(C/) = 0, since at r — > one has 



£iT)^£{0)-aRc{g)pT + 2 



£{0) 



p- 



^ + ^im(t;) 

24 2 ^ ' 



{prf 



(188) 



Formula (188) is exact in the case of 7 = 1. 



If 7 > 1, then one should replace each function sinh(aa;)/a in (185) by its trigonometrical counterpart 
sin(aa::)/a: see Eq. (|5^). In this case the total energy oscillates in time with the period TT/{pd), returning to 
the initial value at the end of each period. For a large detuning 7^1 the amplitude of oscillations decreases 
as 7^^ if Keg ^ and as 7^^ otherwise. 

Note that the total "vacuum" and "nonvacuum" energies increase exponentially with time, if 7 < 1 
and T > 1, whereas the total number of photons increases only as and r, respectively, under the same 
conditions. The origin of such a great difference in the behaviours of the total energy and the total number 



of photons becomes clear, if one looks at the asymptotical formulae (102)-(104). They show that the rate of 
photon generation in the mth completely excited mode decreases approximately as Mm ~ 1/m (excepting 
the modes whose numbers are multiples of p) , so the stationary rate of the energy generation £„i = mMm 
asymptotically almost does not depend on m. In turn, the number of the effectively excited modes in crea ses 
in time exponentially. These two factors lead to the exponential growth of the total energy (see also |163| in 
the special case 7 = 0). 
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8 Three-dimensional nondegenerate cavity 
8.1 Empty cavity 

Now let us proceed to the three-dimensional case. For definiteness we choose a rectangular cavity with 
dimensions L^^ Ly, (briefly designated by symbol {L}). If these dimensions do not depend on time, each 
field mode is determined by three integers m,n,l, responsible for the eigenfrequency 



im/L^f + {n/Lyf + {l/L,f] ^'^ , (189) 



and by two orthogonal directions of polarization. In order to simplify the exposition and to get rid of extra 
unessential indices, let us consider the case when <C L^; ~ Ly. Then the frequencies with I ^ are much 
greater than those with / = 0. It is clear that the interaction between low- and high-frequency modes in 
the nonstationary case is weak. Consequently, studying the excitation of the lowest modes we may confine 
ourselves to the case of / = 0. Then the only possible polarization of the vector potential is along z-axis, so 
the low-frequency part of the Heisenberg field operator at i < reads 

i,(x, y,t<0)^J2 V'n {x, y\{L}) [&„e-''^"* + 6|,e^""*] . (190) 

n 

The difference from the similar expression ( |l^ ) is that now the suffix n is replaced by its "vector" counterpart 
n = {m,n), and the function ip^ {x,y\{L}) depends on two space coordinates: 

ipn[x,y\{L}) = 2[LxLyL,J sm— — sm-— . 

LiX J-iy 



The coefficients in Eq. (190) are chosen again in correspondence with the standard form of the field Hamil- 
tonian (p^. 

Now let the dimension to depend on time according to the given law L{t). To satisfy the boundary 
conditions 

^^L^O ^ ^^\x=L{t) = ^zly=0 = ^z\y=Ly = ^ 



we write the field operator at t > in the same functional form ( |190| ), but with the time-dependent parameter 

L{ty. 

A,ix, y, t) = 2V¥^ ^„ (x, y\L{t),Ly) Q^{t). (191) 

n 

In the stationary case the oper ators (5n(0 coincide with the (coordinate) quadrature components of the field 
mode operators. Putting (191) into the wave equation 

d'^A,/de - AA, = 

we arrive at the equation looking just as Eq. (p^: 

j j 

+A'(i)EMiQ{"^ (192) 
ji 

Now \{t) = L{t)/L{t), all the indices are "two- vectors" , the frequencies are given by Eq. ( |l89| ) with I = 
and L{t) instead of L^, and the constant numerical coefficients g^j are given by the integrals 



gkj^L dx dy dzVj(r|L 
Jo Jq Jo 



dL 



The explicit form of "two-dimensional" coefficients gkj is more complicated than a simple formula (pj|). 
However, these coefficients remain antisymmetrical: gkj = ~5jk, due to the normalization of functions i/jk, 

{L} 
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and due to zero boundary conditions at x = L. (Moreover, they do not depend on the cavity dimensions.) 

Although we use the same notation as in the ID case, the operators Qn{t) in Eq. ( |l9l[ ) differ from 
their analogs in a similar decomposition (|l7|). Now Qn{t) means the Hermitian operator coinciding with the 
("coordinate") quadrature component of the field mode operator. 

Suppose for simplicity that the wall oscillates at twice the eigenfrequency of some unperturbed mode, 

L{t) =La[l- CL cos(2wmO] , IelI < 1 

and let us look for the solution to Eq. ( |192[ ) in the form 

<3k(i) = ^ki^Lt) exp (-it^kO + VkieLt) exp {iujut) (193) 

(we have omitted "hats" over operators). Contrary to the one-dimensional case, now all the terms on the 
right-hand side of Eq. (192) disappear after averaging over fast oscillations, since the spectrum cjj is not 
equidistant. Indeed, the first and the second sums on the right-hand side do not contain functions Qk due 



to the antisymmetricity of coefficients gkj, whereas the last sum is proportional to A 



Consequently, 



after the multiplication by the proper exponential functions, the right-hand side will consist of the terms 
containing the factors such as exp {i [±Wj ± Wk ± "^LOm] t) with j 7^ k. After averaging all these terms turn 
into zero. (Strictly speaking, the frequency spectrum ( |189D contains the "equidistant subset" , corresponding 
to the indices m, n, I multiplicated by the same integral factors. However, this fact does not change the 
conclusion, because the "coupling constants" t/kj between such modes are equ al to zero.) 

Consequently, in the resonance case the field problem is reduced | 188 , 18£] to that of a one-dimensional 
parametric oscillator with the time dependence of the eigenfrequency in the form 



uj{t) = tJo [1 + cos{2ujot)] 



(194) 



ojQ = LOmn being the unperturbed eigenfrequency of the resonance mode. Here the frequency modulation 
depth e is related to the cavity length modulation depth as follows: 



1 



1 -I- {nLo/mLyY 



-1/2 



We use the notation e in order not to confuse the dimensionless modulation parameters in the one-dimensional 
and three-dimensional cases. 

At this point we may abandon the Heisenberg picture and proceed to the Schrodinger representation. 
Of course, both representations are equivalent, as soon as the field problem has been reduced to studying 
a finite-dimensional quantum system. However, the most of numerous investigations of the time-dependent 
quantum oscillator, since Husimi's paper [285|, were performed in the Schrodinger picture. So it is natural 
to use the known results. According to |279,28E, 286|), all the characteristics of the quantum oscillator are 



determined completely by the complex solution of the classical oscillator equation of motion 



satisfying the normalization condition 



. + uj^{t)u = 0, 



2i. 



(195) 



(196) 



Let us assume that function a;(t) takes the constant value loq at t < and at t > tf > 0. Moreover, it is 
convenient to choose the initial conditions for u-function as follows: 



u{0) ^ l/y/uJ^, u{0) = i^/uJ^. 



(197) 



Then the quantum mechanical average number of photons created from the ground state due to the time 
dependence of the frequency in the interval of time < t < tf is given by the formula 



Looking for the solution of Eq. ( |195| ) in the parametric resonance case (194) in the form 

1 



uit) = 



Tj{t)e 



— iujoti 



(198) 



(199) 
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(the opposite signs in the arguments of the exponential functions in Eqs. (193) and (199) are due to the 
different representations: the former equation is written in the Heiscnbcrg picture, while the latter ~ in the 
Schrodinger one) and using the method of averaging over fast oscillations, one can easily obtain the first 
order differential equations for the amplitudes (provided that |e| <C 1), 

^ — iujQsr], 7] = —iujQe^,. (200) 

Their solutions satisfying initial conditions (197) (up to the terms of the order of e) read [|l|,|8^-|90) 

£, {t) = cosh{uJo st) , ri{t) ~ —isinh{u;oet). (201) 



Due to Eqs. (198), (199) and ( pOl[ ), the average number of photons (and the total energy in the cavity) 
grows in time exponentially: 

(n) ^ = sinh^ (wqeO ■ (202) 

It is well known that the initial vacuum state of the oscillator is transformed into the squeezed vacuum 
state, if the frequency depends on time (see, e.g., reviews |2|,|79| and numerous references therein). More- 
over, looking at Eq. (202) one can immediately recognize the combination tooit as the so called squeezing 



parameter. Therefore the probability to registrate n photons exhibits typical oscillations: 



[tanh (c^Qg^)]"" (2m)! 
cosh{uJoet) (2™m!)^ 



2m+l 



(203) 



This distribution possesses the photon number variance cr„ = ^ sinh^ (2LjQet) . Similar formulas for the 
amount of photons created in a cavity filled with a medium with a time-dependent dielectric permeability 
(and stationary boundaries) were found in |222]. The quadrature variances change in time as (now t — sujQt) 



U 



-2t 



2t 



(204) 



An unlimited squeezing can be achieved in this case due to the absence of interaction with other modes. 



8.2 Interaction with a probe oscillator inside the cavity 

The situation changes drastically, if the field mode can interact with some detector placed inside the cavity. 
Following |188,189| we demonstrate the effect in the framework of a simplified model, when a harmonic 
oscillator tuned to the frequency of the resonant mode is placed at the point of maximum of the amplitude 
mode function ipmn {x, y\{L}) in the 3D rectangular cavity. 

Assuming the interaction between the oscillator and the field to be described by means of the standard 
minimal coupling term — (e/TOc)pA, wc arrive at the following two-dimensional Hamiltonian governing the 
evolution of the coupled system "field oscillator -I- detector" : 



H=l [P'' + Lo\t)Q^+p^+LO-iq 



2„2 



(205) 



Here P, Q are the quadrature components of the field oscillator, and p, q are those of the probe oscillator. We 
neglect the interaction with nonresonant modes, since it is reasonable to suppose that under the resonance 
conditions their contribution is not essential at e <C 1. 



In general, the dimensionless coupling coefficient k must depend on time, due to the decomposition (191). 
However, since this coefficient is small, its variations of the order of en can be neglected in comparison with 
the relative variation of the eigenfrequency Slo/lu ^ e. So k is assumed to be constant. 

Suppose that the lowest cavity mode is resonant. Then one can evaluate the dimensionless coupling 

1/2 

constant as k ^ (e^/27rmc^L) (here we return to the dimensional variables). The maximum value of 
parameter e is (see the discussion in section |l0| ) Emax ^ Smax^s/^Trc, where Smax ^ 0.01 is the maximal 
possible relative deformation in the material of the wall, and Vs ^ 5 ■ 10^ m/s is the sound velocity inside 

1 /2 

the wall. Then the ratio s/k cannot exceed the value Smax (mv'^L/8ne^) ^ 0.05 for i '-^ 1 cm and m ^ 
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the mass of electron (for these parameters k ^ 2 ■ 10^^). Consequently, one may believe that in the real 
conditions i/K<^ 1. 

In the time-independent case, uj(t) — const = loq, we have two eigenfrequencies 

w±=cjo(1±k) (206) 

(provided that |/t| ^ 1). Let us assume that the wall vibrates exactly at twice the lower frequency w_: 

uj{t) = cjo [1 + 2ecos(2w_t)] . (207) 

Then the lower and upper modes practically do not i nteract in the limit of e <C k. 

The Schrodinger equation with Hamiltonian (205) can be solved in the framework of the gener al th eory 
of multidimensional quantum systems with arbitrary quadratic Hamiltonians, first proposed in | 278 | and 
exposed in detail, e.g., in |27£]. In particular, if both for the field and the probe oscillators were initially in 
their ground states. 



i/j{Q,q,0) = \/-exp 



(208) 



then the wave function of the coupled "field + probe oscillator" system at t > can be written as [189| 
HQ, t) = Jg^^^ cxp (^-tt - i [ait)Q^ + bit)q^ ~ 2cit)qQ] ^ , 



(209) 



with the following coefficients: 

a{t) = l + i tanh /ie^^*"^- - i^e'* [tanh //e~** (l + tanh fi sin <^e'f) 



sm ip\ 



b{t) = 1 - itanh^e"^"^- - i^e'* [tanh^e"'* (l - tanh ^ sin $e"^) + sini^] 



z{t) = tanh/ie 



1 - COS ipe-'^ + i tanh^ ^ sin ^e'^'^'^*) 



Here 



$ — {uj^ + a;_) t = 2LUot, Lp = [oj^ — CiJ_) t = 2LUoKt ~ K$, 

/i = £Wo^, A = 1 — cosh \i COS 

In all the formulas above, the terms of the order of k} were neglected, as well as the terms proportional to 
e (excepting, of course, the arguments of the hyperbolic functions). Evidently, $ 3> 3> /i. Hereafter we 
confine ourselves to the most interesting long-time limit case, when 3> 1. Then all the terms proportional 
to K can be neglected, so one can write 



a(i)=:l + ?;x, 6(t) = l-ix, c(t) = x, x==tanh^e 



Eq. (209) shows that the coupled system turns out in a two-mode squeezed state at i > 0. The properties 
of this state, as well as of any Gaussian state are determined completely by its covariance matrix 



M : 



\\M 



a/3 I 



Mal3 



1 



{ZaZi3 + ZpZa), 



where the 4-dimensional (in the present case) vector z is defined as follows: z — (tt, x) = (P, p, Q, q) (evidently, 
(z) = in the case under study). Using the general formulas for multidimensional Gaussian states given 
in [^79| , we have obtained |18S| the following explicit expressions for the two-dimensional blocks of matrix 
M in the long-time limit /x ^ 1: 



1 2 
Mtttt = 2 cosh /i 



1 2 

y^xx = 2 cosh ^ 



1 + tanh /i sin <f> — tanh /i cos (j) 

— tanh /i cos 4> 1 — tanh ^ sin < 

1 — tanh sin (p tanh /i cos (p 

tanh cos 4> 1 -f tanh ji sin (; 
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■ sinh(2/x) 



— COS(/) 

tanh /i — sin ( 



tanh /i — sin ( 
cos (h 



where — 2ip-, and M means the transposed matrix. Consequently, there exists a strong correlation between 
the field and probe oscillators in the long-time limit. For instance, the correlation coefficient between the 
quadrature components reads 



(qQ) 



sinh fi cos < 



(If (j) K, 7r/2, this coefficient, as well as other analogous elements of the covariance matrix, does not turn 
exactly into zero; in such a special case VqQ ~ k, due to neglected terms of the order of k.) 

It is clear that the density matrix of the probe oscillator (which is obtained from the density matrix of 
the total system p {Q, q; Q\ q') = ijj {Q, q) ip* (Q' , q') by putting Q = Q' and integrating over Q) has also the 
Gaussian form. Its properties are determined completely by the reduced covariance matrix (accidentally, it 
coincides with M-^x when luq = 1) 



1 



■ cosh fi 



1 — tanh /i sin ( 
tanh /i cos 



tanh fi cos 
1 + tanh fi sin ( 



(210) 



A similar matrix for the field oscillator can be obtained from Eq. ( 210 ) by means of changing the sign 
of parameter p.. As was shown in [280|, the photon statistics in Gaussian one-mode states is determined 
completely by two invariants of the covariance matrix. 



d = detM, 



T = Tr M. 



(211) 



Evidently, parameter T is twice the energy of quantum fluctuations. The parameter d characterizes the 
degree of purity of the quantum state, due to the relation [277| 



2Vd' 



(212) 



where p is the statistical operator of the system. The degree of squeezing, defined as the minimal possible 
value of the variance of some quadrature component, normalized by its vacuum value (which is equal to 
h/{2muj) for an oscillator with mass m and frequency w), is determined jointly by both parameters, T and 
d, according to the relation [p7l [cf. Eq. (|10C|)] 



Both subsystems have identical invariants: 



2(q2) = T-VT^ 



T = 4d = cosh^ /X, 



Ad. 



(213) 



so for /X 3> 1 they appear in highly mixed quantum states, with rather moderate degree of squeezing, which 
tends asymptotically to 50% (cf. the one-dimensional case described in subsection 6.3): 



s = e cosh/_i = 2 ("'^ 



The average number of quanta in each subsystem equals 



1 



1 



(T — 1) = — sinh /X, 



^ ' 2' '2 

i.e., twice less than in the case of an empty cavity. The variance of the number of quanta (photons) equals 



1 



(2T2 - 4d - 1) 



1 . 9 

— sinh" p cosh(2/i) . 
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Mandel's parameter (defined by Eq. (155)) turns out much greater than unity for /i ^ 1, indicating that 
the photon statistics is highly super-Poissonian: 

Q«sinh2(/x). 

The photon distr ibut ion function can be expressed in terms of the Legendre polynomials, according to 
the general formula (151): 

n,, = , .. Pn{~iz), (214) 



where 



3 cosh ^ 
sinh/i 



3 cosh jjL 



Actually the right-hand side of Eq. (214) is a polynomial of degree n with respect to the variable z^, due to 
the recurrence relation 

nVn = 2^ [(2n - l)Vn-i + {n- l)Vn-2] ■ 

li ^ ^ 1, then w 1/3. The behavior of the distribution function (214) was shown in |18£]. Since the 
argument of the Legendre polynomial is pure imaginary, Vn has no oscillations, in contradistinction to the 
vacuum squeezed state. 



8.3 Interaction with a two-level detector 



Another model of the detector, which has only two energy levels, was considered in [188|. The most sig- 
nificant features can be described, in the rotating wave approximation, in the framework of the following 
generalization of the Jaynes-Cummings Hamiltonian: 



1 £ 
H — a^a + —fldz + n (acr+ + a^o"-) + — sm{ujiut) 



(215) 



The eigenfrequency of the unperturbed mode is assumed luq = 1, il is the energy level difference of the 
detector, k and e k have the same meaning as above; a,a^ and CT+,cr_,(T2 are the standard photon and 
spin operators. The wave function of the system "field -|- detector" can be written as 



oo 

m = T.{^n-Ht)\n,-)+ci+\t)\n,+) 



with a clear meaning of the symbols. If e = 0, then the known solution of the JC- model reads [291| 

cVit) = c, (216) 

(217) 
(218) 



cl+i (<) = a„ cos exp - 6„ sin exp [-itE^ ] , 

ci^'' (t) = On sin dn exp [— iii?^] + 6„ cos i9„ exp [—itE"] , 



where rt = 0, 1, 2, . 



E^ = n + -± A„, A„ = 



1/2 



tan 'dn 



2\r. 



1 + 9. 



2A„ + 1 - f7 



1/2 



We suppose that initially the system was in the ground state with the only nonzero coefficient Cg (0) = 1, and 
that the frequency of wall's vibrations is close to twice the frequency of the unperturbed mode: w^, ~ 2 — v. 
Looking for the solution at e ^ in the same form ( ^16| )-(218), but with time-dependent coefficients, and 
neglecting the rapidly oscillating terms containing exp(2ji), we get the following equation for the coefficient 
c(t): 

\/2 

c = — — e {bi sin??i exp [it (Ai — v)] — ai cosi?! exp [—it (Ai -f i/)]} . 
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1/2 

Assuming v = Xi = — 51)^ + 2k^] and neglecting the terms oscillating with the frequencies of the 
order of k, one can check that the infinite system of equations for a„ and 6„ is reduced to the following two 
equations: 

c = ^-esmui Oi, oi — — ^£sm?7ic. 
Consequently, in the resonance case we have only three nonzero amplitudes: 



Cq ■* = cos(at), C2 •* = sinz?i sin(Q!i) exp 



-itEl 



(-) 



.(+) 



cos -di sin(Q!i) exp —itEl 



A-) 



where a = sin i?i£V2/4. No more than two photons can be created, and the probability of finding the 
detector in an excited state is always less than 1/2. All the probabilities, in contradistinction to the 
first example, are periodically oscillating functions of time: 

Vi = 7'^+' = cos^ sin^(a<), 1^2 = sin^ ?9i sin^(ai). 

It is interesting that the upper level of the detector never can be populated with 100% probability, since 
i3„ > for all values of parameters. For = 1, z9„ = 7r/4, and Vi = V2 = V'^^^ = ^siv?{£t/A). Large 

detuning, 1 — 17 3> k, results in increasing P^axi since 0. However, in such a case a — > 0, as well, and 

the applicability of JC-model to the description of the interaction between the detector and field becomes 



questionable. Recently, a more general model was considered in |238|. 



9 Influence of damping 



The complete theory of the field quantization in media with moving nonideal boundaries is not available 
at present. The field quantization in spatially inhomogeneous but nonabsorbing dielectrics was studied 
in 1 292, 293| , and the same problem for nonabsorbing media with time dependent parameters was considered 
in [217-219, 222|. The case of absorbing media was analysed in [294-297|. The theory of the field quantization 



in leaky cavities was developed in [298-B00|. However, in all those studies the boundaries were fixed. Due 
to the complexity of the problem, only a few simplest models have been considered up to now in the case of 
moving walls. 

For example, one can try, as the first step, to neglect coupling between different field modes inside the 
cavity. Such an approximation can be justified, e.g. for an adiabatic motion of the cavity walls, when the 



characteristic mechanical frequency. 



is many orders of magnitude less than the electromagnetic field 



eigenfrequency, u>e- However, no new photons can be created under the condition Um ^ uJe (and the photon 
number distribution cannot be changed, as well), since the photon number operator is the adiabatic invariant 
in this case. 

Fortunately, as was shown in the preceding section, the interaction between different field modes can 
be neglected also in the case of a three-dimensional cavity with a nonequidistant spectrum of the field 
eigenfrequencies, under the parametric resonance condition ujm ~ 20;^. In such a case, one can infere some 
quantitative information on the behavior of the field in the cavity, studying the problem of the parametrically 
excited oscillator with damping. 

The influence of an environment on the parametric amplification was considered in detail, e.g., in |301[ , 
where an explicit coupling with a heat bath consisting of harmonic oscillators was introduced. More ge nera l 
models of the environment were studied, e.g., in the framework o f the inf l uence fu nctional approach |302| , 
mainly in connection with cosmological problems. It was shown in [299,30C, t303 ,304l that the influence of the 
"modes of the universe" outside the cavity with fixed mirrors can be described effectively in the framework 
of the Heisenberg-Langevin equation of motion for the photon annihilation and creation operators a,a^^. 
An equivalent description in the Schrodinger picture is achieved in the framework of the "standard master 
equation" |305-310|, whose simplest form reads 



a)ap - 



pa)a\ 



(219) 
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where p is the statistical operator of the distinguished field mode, H is the Hamiltonian, and the damping 
coefficient 7 absorbs all the details of the loss mechanism: the transmissivity of the mirror, the coupling 
between the field and the atoms inside the wal l , etc ., so that it is proportional to the reciprocal of the 
dissipation time scale. It was assumed in [256, 257 that Eq. (219) can be used as well in the case of 



the moving mirrors , with the same value of the damping coefhcient as in the case of the fixed boundaries. 
Following the same line as in i25^ , p57| , we also assume that the time evolution of the mixed quantum state of 
the resonance field mode is governed by a linear master equation (although one cannot exclude a possibility 
that such an approach is oversimplified: see, e.g., 311 , 312t ). However, instead of using the operator equation 
like ( |2I9| ), we consider the most general linear equation of the Fokker-Planck type for the Wigner function 
W(q,p,t) |277, 30§| , 3I3| ] [q^p are the quadrature components of the field mode). 



dW 



|(M-P]I4^) + | 



D. 



{t)q]W) 



'qq 



dq^ 



dp^ 



qp 



dqdp' 
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The coefficients 7^ and Dij — Dji depend on the concrete form of the microscopic interaction between the 
system involved and an environment |305-30S, 3f4|. For example, the simplest models of the damped optical 
oscillator with a constant frequency ojo yield the following set of coefficients 30^ , 315 1: 



7p = 7? = 7i"s(wo)l5(^o)l^ 



Dpp = LolDqq = 7gfeg(^^o), 



Dpq = 0, 



(221) 
(222) 



where s{uo) is the density of states of the reservoir, g{(jJo) is the function describing the intensity of couphng 
between the distinguished oscillator and the reservoir degrees of freedom, and Eeq{u}o) is the equilibrium 
energy of the oscillator with frequency ojq at temperature T, 



^ea (^0) = —flUJn COth ( — — %- 
eqy 0/ 2 " \2kBT 



(223) 



Choosing different couplings between the oscil lator und er study and the reservoir, one can obtai n variou s 
other sets of the drift and diffusion coefficients |30S , 314 1 , but all of them must obey the constraint 1 277 , 3I£ - 



D D — 



> 



fl'^hp+lq) /16, 



(224) 



which guarantees an absence of nonphysical solutions violating the uncertainty relations and corresponding to 
nonpositively definite density matrices. Besides, the coefficients Dpp and Dqq must be positive. However, it 
will be shown in this section, that in the case of a weak damping, the evolution depends on two combinations 
of the damping and diffusion coefficients only. 



7 = ^ (7p + 7?) 



f * — — {Dpp + LolDqq) 



(225) 



Due to the fluctuation-dissipation theorem, the diffusion coefficients are proportional to the damping coef- 
ficients, therefore does not depend on 7, at least up to small corrections of an order of 7^. The physical 
meaning of the parameters 7 and is elucidated below: 27 is the reciprocal energy relaxation time of the 
cavity due to all possible mechanisms (a real dissipation in the walls and the leakage through the boundaries), 
i.e. 27 = loq/Q, Q being the cavity quality factor, whereas f» = £eq i^o)- 



Solutions to Eq. (E20[) with a constant frequency and different sets of constant diffusion coefficients were 



obtained by many authors; they were analysed in detail in |277, 308 1, where other references can be found. 
Since Eq. ( |22C| ) looks like a two-dimensional Schrodinger equation with a quadratic (although nonhermitian) 
Hamiltonian, its propagator can be calculated in the framework of the method of quantum integrals of 
motion p78[. The explicit form of this propagator in the generic case of time dependent coefficients can be 
found in | |27^ , p20t . Here we confine ourselves to calculating the second order statistical moments and the 
energy (the number of photons) of the field oscillator. 
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9.1 Evolution of the energy and the second order moments 



The time dependence of the energy and the second order statistical moments (variances) of the field mode 
quadrature components, aab = ^{ab + ha) — {a){b), is governed by the equations following from Eq. (122C 



(J. 



(Tnn = app — U}~[t)aqq — 'Zjann + 2D. 
O'pp — SU-' (t)(Tpq 2,^p(Jpp -r jLi^pp 

We assume that the oscillator eigenfrequency depends on time as 



= 2apq — 2^qaqq + 2Dqq, 
= Gpp — U?{t)(Jqq — 2japq t ^-L^pq, 

2u\^ ^ I on 



UJ 



it)^u;o[l + 2esm{nt)], n = 2{ujQ + S), < wq, |e| < 1, 



(226) 
(227) 
(228) 



(229) 



where loq is the unperturbed field eigenfrequency, and is the frequency of the wall vibrations. Then one 
could suppose that the damping and diffusion coefficients must depend on time, as well. We argue, however, 
that in the case under study the coef ficient s and Dab can be considered as time independent. For example, 
let us look at the expressions (221) and ( |222| ). In the case of a vibrating cavity, the time variable could 
enter the coefficients ja and Dab through the coupling function g{uJo), which can depend on the variable 
length of the cavity L{t) = Lq [1 + ^Lesin(Sli)], being a numerical coefficient. Then one could expect a 
similar time dependence of the coefficients of the Fokker-Planck equations, j{t) = 70 [1 + £,-y£ sin(Sli)] , D{t) = 
Dq [1 + S^ds sin{Qt)], and being some other numerical coefficients. One should remember, however, 
that the modulation parameter e is very small under the realistic conditions: its absolute value cannot exceed 
The set of equations (226)-(228) contains three small dimensionless parameters, e. 



10~** |188| , |189| . The set of equations (226)-(228) contains three small dimensionless parameters, e, S/uq, 
and 7o/wo, and we are interested in the weak damping case, when these parameters are of the same order of 
magnitude. Under this condition, the time dependent parts of the coefficients 7a and Dab are proportional 
to the products e^q/ujq ~ 0{£^) ~ 10~^^, so it seems reasonable to neglect these extremely small terms. 

A significant time dependence of the damping and diffusion coefficients could arise in the specific case 
of an unstable reservoir, provided that the reservoir oscillators having the frequencies close to uoq could be 
also excited due to some resonance processes between the vibrating surface of the wall and the reservoir. In 
such a case, the energy of the resonant oscillator s wo uld increase in time, resulting in increasing values of 
the damping and diffusion coefficients (see Eq. (|222| )). As a consequence, we would obtain an additional 
amplification of the energy of the field mode due to the interaction with the reservoir. However, such a 
model seems unrealistic, since it implies that some distinguished degrees of freedom of the reservoir are 
practically isolated from the rest of the reservoir. This conjecture contradicts the usual concept of the 
reservoir consisting of a great number of strongly interacting particles, so that the state of the reservoir is 
not sensitive to small external perturbations. Therefore, we assume that the only time dependent coefficient 
in Eqs. (p26|)-(p28|) is w(t). 

The set of equations ( |226| )-( |22^ ) is equivalent to the third order differential equation for the variance 

(Jqq = a 



. da 



IF ^ ^^IF ^ '^'^^^^^'i ^ ^''^^'^ ^ 27tJo)'7 = 875* 
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where the coefficients 7 and 8^, are given by Eq. (225) (we neglect the terms of the second order with 
respect to 7, Dab^ and e). To find an approximate explicit solution to Eq. (230) with function Lo[t) given 
by Eq. ( |229|), we use the well known method of slowly varying amplitudes, which was exposed, e.g. in 
textbooks |266- ^68 | and applied to the quantum parametric oscillator in [ 289| , 301|. Following this method, 
we write 

aqq{t) = A{t) + B{t) cos(rJi) + C(t) sin(m). (231) 



so that the function (231) with A,B,C = const is an exact solution to Eq. (230) for e — ^ ~ 5 ^ Q. 
Supposing that the dimensionless small parameters ^ /lO{) and have the same orders of magnitude as 

the small parameter e, i.e. 7 = 7£u;o, 5 = Setoo, 7, (5 ^ ^(l)j we assume that the amplitude coefficients 
A,B,C are functions of the "slow time" t — it, so that the time derivatives A/ dt^ ,d'' B / dt^ ,d''c/ dt^ 
are proportional to [k = 1,2,3). Then we put the function (231) into Eq. (230) and neglect the terms 
proportional to and e^. Besides, we perform averaging over fast oscillations with the frequency il, in 
order to eliminate higher harmonics with frequences mfi, m — 2,3, . . ., whose amplitudes are proportional 
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to e™"^ 1 267, 268 1). Finally, we arrive at the following set of equations for the slowly varying amplitudes 
(here the overdot means the derivative with respect to the real time i, and the new parameter k = euo has 
the dimensionality of frequency, like 7 and 6): 



A = -2jA + 2kB + 2'y£^, 
B = 2nA - 2-fB - 25C, 
C = 26B~ 27C. 



(232) 
(233) 

(234) 



This system can be easily solved for an arbitrary time dependent function 7(t), since the substitution 
A{t) = A(t) exp [—2 / j{T)dT\ removes the function j{t) from the homogeneous parts of Eqs. ( |232| )-(234). 
However, we consider the case of constant coefficients only, due to the physical reasons discussed above. 

To simplify the formulas, we assume hereafter h — loq — I; thence Q = 2 in the amplitude coefficients. 
Neglecting small terms of the order of 7/^^, n/fl, and in the amplitude coefficients, one can express the 



variances app and Upg as follows 



app{t) = A{t) - B{t) cos(m) - Cit) sin(f7t), 
[t) = C(t) cos(m) - B(t) s\n{VLt). 



(235) 
(236) 



Then the initial conditions for the set ( |232| )-( ^34| ) read 

A{0) = [a,,(0) + app(O)] /2, 5(0) = [(t,,(0) - (Tpp(O)] /2, C(0) = ap,{Q). 

Evidently, A{t) coincides with the mean energy to within small corrections of the order of e, 7, and 5: 

E{t) = {<Tpp + aqq) /2^ A{t). 

In order to elucidate the meaning of the coefhcients B and C, consider the determinant of the invariance 
matrix (cf. Eq. (|2Tl| )) 

d = appaqq ~ al > h^/A. (237) 



Here the last inequality holds due to the Schrodinger-Robertson uncertainty relation [271, 27!; ,321, 322 1. The 
meaning of the parameter d as the universal quantum invariant is discussed in J323| . The minimal invariant 
variance (|100|) can be expressed in terms of £ and d as 



u = E - ^/£^ -d = 

and one can easily verify the relations 

d^ A^ - B"^ -C"^, 



E + VF^~d 



The solutions to Eqs. (232)-(234) read 
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= A, - 




-27* 
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= C>^ ^ 
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where v = \J — i^^, and 
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a± 
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-4 [kC(0) - M(0) + EJ] , 
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(239) 
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The meanings of the parameters 7 and become clear, if one considers the special case of the oscillator with 
time independent coefRcients, k ~ 5 = v = Q. Then we see that 27 is the energy relaxation coefRcient, so that 
it can be expressed in terms of the cavity Q-factor by means of the relation 27 ~ loo/Q- The energy of the 
oscillator in this special case tends to as t ^ cxd, and this value can be identified with the thermodynamic 
equilibrium oscillator energy £eq given by Eq. (223) (up to corrections of the order of (7/1^0)^ |277, 3081. 

The sign of the difference v — ^ determines the regions of stable and unstable solutions of Eq. (23C) in 
the space of parameters k, (5, 7 (for small values of these parameters). The stable (limited in time) solutions 
exist for large values of the damping or detuning coefficients, v < ^, i.e. 



2 I r2 ^ 2 

7 + d > K 



(245) 



In this case, the final state of the oscillator does not depend on the initial conditions. The asymptotical 
values of the energy, the c?-factor, and the minimal invariant variance read 



£(00) 



7' + ^' 

72 -1.2 --9' 



d{oo) 



1^ + 5^ 



,2 e?7 



u[oo) 



--£e 



(246) 



At zero temperature, when n tends to the threshold, the minimal variance u goes to the value 1/4, which is 
twice less than in the coherent state. 

Above the threshold, i.e. in the instability region v > ^, the energy (or the number of photons) increases 
exponentially for {v — "f)t ^ 1, 

£{t) = a+ne^^"-^^', (247) 



and it depends on the initial conditions through the coefficient a+. Using Eq. ( 244 ) and taking into account 
the uncertainty relation d > 1/4 (for ?i = 1), one can verify that the coefficient a+ is bounded from below 
by a positive value, 

KjSeq 



2l/2(zy-7)' 



Consider a special case of initial thermal equilibrium state. Then 



k£. 



eq 



^0 — 



0, d{Q)^£%, 



thus d-factor ( ^39|) depends on time as 



d{t) 



E. 



eq 



viv'^ ~ 72) 



(248) 



(249) 



If 7 > and [v — 3> 1, then due to Eq. (238), the minimal variance tends asymptotically to a constant 
value 

^ ^ (250) 



-£e 



"f + v 

Consequently, a large squeezing can be achieved even for a high temperature initial state, if 3> 7. If 7 = 0, 
then d does not depend on time, d = £g^, and for vt ^ \ the minimal variance goes asymptotically to zero as 
u w £eq{v I exp{~'^vt)- Ouc should remember, nonetheless, that the solutions indicating the exponential 
growth of the energy are justified until t t2 ^ {uJoe^)~^, since for larger times the neglected second order 
terms in Eqs. (23C), ( |232| )-(234) could become important. However, the time t2 is very large under the 
realistic conditions. 



10 Discussion 

We have demonstrated a significant progress in our understanding and quantitative description of quantum 
processes in cavities with moving boundaries, achieved 30 years after the pioneer paper by Moore [ [76| and 
almost 80 years after the first papers on classical electrodynamics in such cavities by Nicolai and Havelock 
l^-H). We have shown that quanta of electromagnetic field can be created from vacuum in a cavity with 
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vibrating walls under the resonance condition (and cannot be generated for nonresonance nonrelativistic laws 
of boundary motion, in particular, in the case of a large detuning from the resonance), and the quantum 
state of field exhibits the "nonclassical" properties. 

The possibility of observing the effect depends crucially on the achievable values of the wall displacement 
amplitude. For the cavity dimensions of the order of l-;-100 centimeters, the resonance frequency loq/-k varies 
from 30 GHz to 300 MHz. It is difficult to imagine that the wall could be forced to oscillate as a whole at such 
a high frequency. Rather, one could think on the oscillations of the surface of the cavity wall. In such a case 
one has to find a way of exciting a sufhciently strong standing acoustic wave at frequency ojw — 2^0 inside 
the wall. The amplitude a of this wave (coinciding with the amplitude of oscillations of the free surface) is 
connected with the relative deformation amplitude 5 inside the wall as J = ujyi;a/vs, where Vs is the sound 
velocity. Since the usual materials cannot bear the deformations exceeding the value S^ax ~ 10^^ [ |324| , the 
maximal possible velocity of the boundary appears Vj^ax ^ ^max^s ^ 50 m/s (independent on the frequency). 
Thus the maximal dimensionless displacement e = a/Lg is e„iax ~ {vs/2Trc)6max ~ 3 • 10~* for the lowest 
mode with the frequency o^o ^ ctt/Lq. It also does not depend on the frequency. Consequently, the maximal 
rate of photon generation in the principal mode of a ID cavity can be estimated as 



, 2 '^ma.TT -6- 10- V/27r. (251) 

J max ^ C 2tT 

It is proportional to the frequency. For uji/2tt — 10 GHz (corresponding to a distance between the plates of 
the order of several centimeters) we get 600 photons/sec. 

This number can be significantly increased in a 3D cavity, due to the exponential law ( ^02[) . For the same 
frequency loqI^'k = 10 GHz, the maximal value of parameter fi — ^toot equals fimax ~ 600t, time t being 
expressed in seconds. Even if the amplitude of the vibrations were 100 times less than the maximal possible 
value, in i = 1 s one could get about sinh^(6) ~ 4 • 10** photons in an empty cavity. Obviously, the concrete 
shape of a 3D cavity is not important. The significant requirements are: i) the nondegenerate character of 
the eigenfrequency spectrum, and ii) the condition of the parametric resonance between the oscillating wall 
and some electromagnetic mode. The total energy of photons created in the ID cavity is approximately the 
same as in the three-dimensional case. The difference is that in the ID case this energy is spread over many 
interacting modes, resulting in moderate numbers of quanta in each mode. The rate of photon generation in 
the m-th (odd) mode of the ID cavity is approximately m times less than in the principal mode with m = 1 
(in the asymptotical regime). 

To create the above-mentioned 600 or 4-10"' photons, one should vibrate the wall for not less than 1 second. 
The necessary Q-factor of the cavity must be Q ^ 3 • lO^*'. This value was achieved in experiments already 



several years ago |325]. An unsolved problem is how to excite the high-frequency surface vibrations with a 
sufficiently large amplitude. One could think, for instance, on using some kind of piezoeffect. This method 
was successfully applied in early experiments devoted to solving the mode-locking and pulse production 
problems in lasers with the aid of vibrating mirrors. The displacements of the mirror from 0.1 /im to 0.7 /im 



at the frequency 500 KHz were achieved in |27, 28[. In |Q the resonance vibrations of the mirror in a laser 
with the length 250 cm (i.e., at a frequency about 100 MHz) were excited with the aid of a quarz transducer. 
However, for our purposes the frequency 100 MHz is too small, since the parameter fi becomes 100 times 
less, comparing with the estimations given ab ove (rem e mbe r that Smax does not depend on the frequency). 
Fortunately, the results of recent studies [192, 21C| -212] show that the influence of temperature is not 



so significant as it could appear at first glance. Moreover, in certain situations the initial temperature 



fluctuations could be used to amplify the effect |211|. However, the resonance requirements are rather hard: 
for the values of the frequency uj ~ 10^*^ Hz and the maximal possible modulation parameter e ~ 10"* the 
admissible detuning should not exceed 100 Hz during the whole time interval ^ 1 s, necessary to accumulate 
the resonance effect. 

One of the reasons of the studies on the dynamical Casimir effect for the last few years was Schwinger's 



hypothesis |152-157| that this effect could explain the sonoluminescence phenomenon, i.e., the emission of 
bright short pulses of the visible light from the gas bubbles in the water, when the bubbles pulsate due to the 
pressure oscillations in a strong standing acoustic wave. (For the reviews and numerous references related to 



this effect see, e.g., 121, 326- 328| |.) There are several publications, e.g., |329-331], whose authors considered 



the models giving tremendous numbers of photons which could be produced even in the visible range due to 
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the fast motion of the boundaries. However, the analysis of these models shows that they are based on such 
laws of motion of the boundaries which imply the supcrluminal velocities, so they are not realistic. 

Although the results of this chapter, being obtained in the framework of simplified one-dimensional and 
three-dimensional models, cannot be applied directly to the analysis of the sonoluminescence problem, they 
are not in favor of Schwinger's hypothesis. The main difficulty is connected with quite different time scales 
of the phenomena. The accumulation of the 'dynamical Casimir energy' is a very slow process, which needs 
a great number of wall oscillations, whereas the sonoluminescence pulses (containing up to 10'' photons) 
have the duration of the order of picoseconds. Moreover, the wall oscillations must be in extremely fine 
tuned resonance with the field eigenfrequencies, since the detuning S > e completely destroys the energy 



growth 1 191]. In particular, if the frequency of the wall oscillations LUwaii is much less than the minimal 
field eigenfrequcncy lui, then the field variation is adiabatic, and the mean number of created photons is 
proportional to [189| {uj^aii/^i)'^ *C 1- These features survive in the three-dimensional model considered in 



this chapter, too. Therefore, it is difficult to believe that very specific conditions of the parametric resonance 
described above could arise naturally in the sonoluminescence case. For other discussions of the problem 



see, e.g., [332| and the contributions in Ref. |333] 



Actually, the main obstacle to produce the 'Casimir light' is the very low ratio of the wall velocity to 
the speed of light in possible laboratory experiments. If the velocity of the boundary were of the order of 
c, then a sufficient number of photons could be created from vacuum practically for any law of motion. For 
the nonrelativistic velocities the only possibility is to accumulate the effect gradually under the resonance 
conditions. Nonetheless, perhaps, the experimental situation could be improved in t he ca se of using some 



kinds of 'effective mirrors', such as, e.g., the layers made of the electron-hole plasma |209|, or some others. 
Therefore, we cannot exclude a possibility that in a not very remote future one could assist a show of 
"quantum magics," when some "quantum magician" takes an empty box, then shakes it well, opens — and 
an astonished audience would see a great number of photons which have appeared "from nothing" due to 
the Nonstationary Casimir Effect. 
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